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ABSTRACT 


Four exponential growth functions: Richards, Brody, 
Bertalanffy and Logistic were fitted to unadjusted data from 
Six breed populations and to data adjusted for year and age 
of dam effects for two populations. Growth curves) and 
parameters characteristic of each function were derived for 


each set of data. 


Data for each breed population consisted of weight/age 
records over ae period of 10 years from 1962-1971 obtained 
from the University of Alberta experimental breeding ranch 


in Kinsella. 


The fit of the functions derived for each set of data 
were examined at birth, weaning, yearling, 20 months, 32 
months, 44 months, 56 months (adult) and over all ages using 
R@ values, mean prediction error estimates and residual 
variances. The degree of fit of the models in the two major 
populations the Hybrid (HY) and Hereford (HE) breed groups 


were studied in detail. 


The Richards and Brody models described growth well 
over the entire trajectory as the R@ values were high while 
mean prediction error estimates and residual variances were 
low. Adjusting the data for year and age of dam reduced the 
residual variances in all models and mean prediction error 


estimates in all but the Richards model. 
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The Richards model provided an accurate prediction of 
birth weight in all of the data set while the Brody model 
provided a variable fit. The Bertalanffy and especially the 
Logistic function were poor in predicting birth weight in 
all data sets. Weaning weight was almost always predicted 
well by the Brody and Richards models. During the 
intermediate ages from yearling to 44 months all functions 
predicted weight with relatively small errors with the 
exception of weight at 32 months of age at which time no 
function could adjust for a temporary weight loss due to 
winter stress characteristic at that age group. The 
functions repeatedly underestimated adult weight and only 
the Richards and Brody functions predicted adult weights 


With some accuracy. 


Adjusting the data for year and age of dam 
Significantly improved the fit of the Brody, Bertalanffy and 
Logistic functions at birth and the Richards at weaning in 
both HY and HE breed groups, while the fit of the Brody was 
improved at weaning in the HY breed group. fhe improved fit 
of the functions by adjusting the data were confined to ages 
prior to yearling and responses were not identical among 


breeds. 


The Richards and Brody functions were consistent in 


predicting small and similar mean prediction error estimates 


» golstvnny eee as hebiwesy ‘absieag 
tehon syhowe ets ae aae ps shee ie a e 
Badia aee 2 at Re elie tend & 


oy pibitergio pita | (wea | o 
ez hipete se aiih pxtoaabo idk fone “ogee 7 . 1 bts yong) 
hiedwt fea each Be geanis io Waniar pakaeen me 
aid ) oe kad ageboe Spee ‘pas via Poser 
eto Se Eis aiyoon ab os paLlanyt: for aatipl e ere 
oat Staw . etoris tame yiiebrstes cod pee 
wi rere etiin, 96" 3 ap 5,20: ataage ae ap Saint 
et bb Pons duper yatzodnes © | 
eae agp eay bape $d, #5) 
yito tne ripen He ieiA eet diame 
ghee taba eegedberg  eaordount 


ve 


fap) 20. eRe DNR ay ‘702. soeb Bae jag 


se ill ons, 3%: $h3 eit phidy ‘ yaa aire boot! ee spa co 
+h bevowqnl oplt Oty bowid: {0 at} ab, ‘padhiew, ea 
eee 9) arp ied ater, Wied ae qithysnthie bi sae te 
“poe Hace soi oe aa one eateeey ade 


ay 
aie i 


pate fdans sa0 saotiasy Be tf axoouk bas ebreday i o as ie 
fas sai Said teat. aoldeae Soe ELaae es 


_ - 


“i. | 


v 


at each age, compared to the Bertalanffy and Logistic 
functions. The overall mean prediction error estimates 
based on the Richards function were small and uniform 
between breeds thus demonstrating an overall consistency of 
fit and was therefore considered as being ideal for 


comparing estimates of fitted parameters. 


The asymptotes (A) fitted by the Richards and Brody 
functions were always higher than those fitted by the 
Bertalanffy and Logistic. The Bertalanffy and Logistic 
functions were characterised as converging rapidly at the 
asymptote as predicted aduit weights and asymptotes were 
Similar, Adjusting data reduced the asymptotes in both 
breed groups, the effects being marked in the fitted 


asymptotes of the Richards function applied to the Hybrids. 


Based on the estimates of the maturing rate parameter 
(k) the Herefords matured faster than the Hybrids. The 
maturing rates of the other breed groups appeared to be 
proportional to the amount of British beef blood in the 


populations. 


Correlations between the parameters (A) and (kK) were 
always negative. Selection for increased growth rate at 
early ages (6 months and 12 months) should result in 
increased maturing rates. Early maturing animals with rapid 


growth rates prior to 1 year of age grew to small mature 
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weights, and late maturing animals with relatively higher 
growth rates at 18 months grew to larger mature weights. 

Animals heavier at inflection showed higher growth rates at 
early ages and lower growth rates at 18 months. Selection 
for higher absolute growth rates at some age would improve 
gain more at immediately adjacent ages as a correlated 


response and less when the ages were widely separated. 


The four growth models could not accurately predict 
absolute growth rate over 4 periods as the observed gains 
were of a fluctuating type and the functions were of an 


exponential decay type. 


The Richards and the simpier Brody models’ were 
particularly good in describing growth in all populations 
and are recommended as being useful and a knowledge of their 
predictive and descriptive properties promotes a better 


understanding of the growth process. 
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1. INTRODUCTION 


Growth in animals is in general the increase in body 
masS with time, and is influenced by a number of factors 
such as the genetic potential of the animal, its nutrition 


and environment. 


An animal basically grows in two ways, either by cell 
proliferation in the tissues (hyperplasia) or by ceil 
enlargement (hypertrophy). The former process assumes much 
importance during the pre-natal growth phase while the 
latter process along with intercellular accretion are 
regarded as major determinants of growth during the post- 


Natal growth phase (Robinson, 1971). 


Growth curves are weight-age curves which are usually 
plotted from birth to maturity on individuals or populations 
and are useful for describing growth patterns characteristic 
of biological types. A gross growth curve can reveal stages 
at which animals make rapid gains in weight as against 
periods when growth is slow or constant, (aS Shown by Brody, 
1945) depending on the slope of the growth curve between two 


points on a time scale. 


Growth curves can also be expressed in other forms such 
as cumulative, absolute or relative and "reflect lifetime 
inter-relationships between an animals' impulse to grow with 


age in all tissues and the environment in which these 
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impulses are expressed" (Fitzhugh, 1976). The idea of 
growth ie. change in biomass with respect to time, was 
given a mathematical description by Brody (1945) in his 
book, "Bioenergetics and Growth." According to Brody a 
cumulative growth curve followed a strict sigmoid shape with 
age and although variations occurred within ages due to the 
effect of the environment, the general shape was 
characteristically sigmoid. Brody divided the weight-age 
curve into three sections, 1) pre-natal, 2) post-natal self 
accelerating, and 3) post-natal self inhibiting phases; the 
latter two being described by two different equations. 

During the post-natal self accelerating phase, growth, 
measured by an increase in body weight, showed an 
exponential relationship with time, Aven chee ae at an 
accelerated rate. During the self inhibiting phase the 
growth process was slower and gains were made at a 
decreasing rate. The point, or period, where the self 
accelerating curve joined the self inhibiting curve was 
recognised by Brody as the point of inflection (POI) and 
considered as the point or stage of sexual maturity 
(puberty). In subsequent studies this point at which the 
velocity of growth changes was computed as being 
approximately the time an animal reaches one third of it's 


mature body weight (Fabens, 1965). 


Following Brody's studies many other equations have 
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been used and developed among which the more important are 
the Gompertz, Logistic, Bertalanffy and the Richards. These 
functions have been applied to study patterns of growth in 
animal forms. Each growth equation has its Own 
characteristics and therefore may give a better fit on 
certain types of data. Most reports that have used a growth 
curve approach concentrate on the overall fit of the 
functions and correlations between fitted parameters, but a 
careful analysis of the fit of a growth curve at various 


ages has not been reported. 


Therefore the present study was designed to evaluate 
the degree of fit of the growth equations - Bertalanffy, 
Brody, Logistic and Richards - to beef cattle growth data 
obtained from breed populations maintained at the University 
of Alberta ranch in Kinsella, with emphasis on their ability 
to describe growth and predict weights during certain pre- 
determined ages, some of which are of economical importance. 
The raw data as well as data adjusted to remove certain 


within age variation were used in the analyses. 


The fitted parameters from each function were compared 
between populations and within functions to arrive at 
estimates of the asymptotes and maturing rates. 
Correlations between parameters both fitted and derived were 
calculated on the functions of best fit to study expected 


responses to selection. The average daily gains based on 
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the predicted weights from each growth function were 
compared with observed average daily gains to determine the 


efficiency of growth curves in predicting this trait. 
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2e LITERATURE REVIEW 


A review of literature on growth curves reveal 
basically two types of equations. The first type of 
equation, which is actually a part whole relationship 
(allometry), was used by Huxley (1932) where (Y) was equal 
to px* in which (b) and (k) were constants and (Y¥) the 
weight of some organ (part) and (X) the body weight (whole). 
Huxley (Cheek, 1968) considered growth as a process of self 
multiplication and the rate of growth of an organism growing 


equally in all parts was at any moment proportional to the 


size of the organisn. 


A second type of growth curve was demonstrated by Brody 
in 1932 (Weinbach, 1941) to be of an exponential nature and 
applied to human growth from birth to puberty. Since the 
publication of Brody's book "Bioenergetics and Growth" in 
1945 many other exponential equations describing weight-age 
relationship were developed, four of which were the 
Bertalanffy (Bertalanffy, 1957), Richards (Richards, 1959), 
Logistic (Nelder, 1961) and the Gompertz (Laird, 1965). 
These functions have been commonly used to describe growth 


of many animal species. 


2.1 Brody function 


The first exponential equation to be developed, given a 
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mathematical description and applied to various animal 
Species was by Brody (1945). Two equations were derived, 1. 
the self accelerating and 2. self inhibiting. The self 
inhibiting equation of the Brody function is a 3 parameter 
model with (A), (B) and («k) as the fitted parameters, 
referred to as growth constants jBrody, 1945). The self 
inhibiting equation is correctly applied post inflection and 
has wide applicability in describing growth of different 


species (Brody, 1945). 


The two equations are: 
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A(1* Be-kt) 
where, | 
Y = weight in relation to time, t. 
a = ¥ intercept 
e = log base e 
B = integration constant established by the values of 
Y andt and adjusts for the situation when 
Y#0;3 
A = asymptote as t > @- 


k = instantaneous relative growth rate (Brody, 1945) 


In equation 1 the growth rate (k) was computed in relation 


to the gains already made hence (k) was positive, and in 
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equation 2 growth rate was computed in relation to the gains 


yet to be made (A+Y) and (k) was negative (Brody, 1945). 


The Bertalanffy function is another 3 parameter 
exponential function used to describe growth (Bertalanffy, 
SESH AD 

Y = A(1 ~ Be-Kt)s 
where, Y, B, A, k and e are as described in 2.1. The 
function is similar to the Brody function and has a fixed 
point of inflection at 8/27 A (Bertalanffy, 1957). Fabens 
(1965) described the Bertalanffy curve as being of a 
decaying exponential type and considered it to have a wide 


application in biology. 


The Richards function was originally an extension of 
the Bertalanffy function and was first applied to describe 


plant growth (Richards, 1959). 


The function is given as follows: 

Y = A(1 - Be-kt)M 
where, Y, B, A, k and e are as described in 2.1 and M = the 
variable shape parameter. Like the other growth functions 


the Richards contains the three parameters (A), (B) and (k) 
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and in addition a fourth parameter (M) which is referred to 
as a variable shape parameter. The value of (M) can vary, 
and determines the point of inflection in relation to the 
asymptotic weight (A) (Richards, 1959). The Richards model 
is a generalized function and the Brody, Bertalanffy and 
Logistic function have been considered as special cases of 


the Richards (Fitzhugh, 1976). 


2-4 Logistic function 


The Logistic function is a 3 parameter exponential 
function and the procedure of fitting it is described by 
Nelder (1961). The function is similar to the Brody 
function and has the following format. 

Y = A(1 + Be-kt)-1 
where, Y, B, A, k and e are as described in 2.1 The Logistic 
function contains the parameters (A), (B) and (k) and its 
point of inflection is fixed at 0.5 A. Its use has been 


demonstrated on mice data by Rutledge et al. (1972). 


2-5 Gompertz function 


SS aoe ae 


Laird (1966) has shown that the post-natal growth for a 
variety of animals such as birds and mammals undergoes an 
exponential decay in specific growth rates with time. The 
Gompertz function has been used to describe growth of mice 


(Laird and Howard, 1967; Laird, Tyler and Barton, 1965; 
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Laird, 1965), and reported to fit better than the Logistic 
function when applied to growth of dairy cattle (Vieira and 


Hoffman, 1977). 


2.6 General 


Growth functions can be fitted to many types of data 
classified as static, cross sectional, longitudinal, mixed 
cross sectional or mixed longitudinal (Cock, 1966; Tanner, 
1951). Data of a longitudinal type provided a relevant 
description of the patterns of growth, as weight for age 
data is recorded throughout an animal's life and a complete 
set of experimental measurements is made available for every 
individual at each stage or agew. In addition to providing 
the data yielded in a static or cross sectional type, it 
provides information on individual variations in growth 
(cock, 1961; Fitzhugh, 1976). This type of data was 
recognised as being most appropriate for describing growth 


in a growth curve analysis (Fitzhugh, 19/76). 


In summary, the basic fitted parameters derived by the 
| parameter models, Brody, Bertalanffy and Logistic 
functions are (A), (B) and (k). The Richards function a 
four parameter model also computes a value for (M) which 


determines the point of inflection. 
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2.7 Application of growth functions 


The Brody function has been used to describe and 
predict weights, and compare growth parameters of many 
species of animals and its wide application in biology has 


been documented (Brody, 1945). 


Brown et al. (1972ab, 1973) fitted the Brody equation 
to female Hereford and Angus growth data and obtained 
estimates of mature weight and maturing rates. The authors 
stated that although the growth curve smoothed out the 
fluctuating data points, the function was adequate in 
describing growth. Their results were not difficult to 
interpret although the genetic and phenotypic correlations 
in relation to the parameter (k) differed considerably 


between breeds. 


The pre- and postrinflection equations of the Brody 
function were used to study growth patterns in Black tailed 
deer (Bandy et al., 1970). They also determined the point 
of inflection in races of Black tailed deer as the pre- 
inflection function followed a straight line up to the point 
of infiection. The authors guestioned the adequacy of the 
Brody function in describing growth rates in absolute terms 
but were satisfied that it was efficient for inter-racial 
Spacer aeons. Lerner (1939) stated that “the pre-inflection 


growth curve although widely used with many excellent 
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results is limited in its value to small areas of a growth 
curve" (Baker, 1944). Brown (1970) described the Brody 
function as one which provided an excellent fit to the data 
points in beef cattle but did not permit an in depth study 


of growth properties. 


Parks (1970ab) developed and used an extension of the 
Brody function: a) relating the biomass of an animal to feed 
consumption, b) cumulative feed consumption to age and c) a 


combination of the first two equations, biomass to age. 


Taylor (1968) used the Brody and Gompertz growth 
equations on data described by Brody (1945) and Laird (1966) 
to demonstrate the proportionality between the time a 
species takes to mature and its mature weight. The time a 
species takes to mature was computed as the 0.2/th power of 
its mature weight; the two growth equations and data giving 
Similar values (Taylor, 1968). Comparing different species 
in relation to their maturing rates Taylor (1968) found 
birds matured faster than mammals, which tends to support 
his statements in an earlier paper relating species size and 


time taken to mature. 


The Richards function plus the Bertalanffy, Logistic 
and Gompertz (Fox, 1970) have been applied to study growth 
patterns of ewes. The Richards function was considered best 


in describing growth although the parameter (M) which 
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determines the point of inflection was negative, thereby 
limiting the study of the inflection parameters. Brown et 
al. (1976) and Brown (1970) found the Richards was better 
than the Brody, Gompertz, Bertalanffy and Logistic functions 


in providing a good fit to beef cattle data. 


The Gompertz function was cited as being adequate in 
describing growth of beef cattle (Cartwright and Joandet, 
1969). The authors stated that the parameters derived had 
biological meaning and could be related to other traits. 
Laird et alj (1966) and Laird and Howard (1967) used the 
Gompertz function to describe growth data for mice and 
analyse strain differences. The authors state that the 


Gompertz function described growth adequately. 


Growth curves of the Logistic and Gompertz type have 
been applied to human height data and a close correspondence 
between observed and predicted values found (Marubini et 
ale, 1971). The Logistic function was applied to study the 
growth among mice selected for low and high body weights 
(Eisen et ale, 1969). The authors reported that the 
Logistic function provided the best fit and description of 
growth among lines of mice although the other equations that 
were tried could have detected differences between lines. 
Rutledge et al. (1972) compared the Logistic and Richards 
functions as applied to mice and concluded that the former 


offered a better fit to the data. A growth function of the 
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Logistic type has been applied to broiler chickens and 
appears to be efficient in describing growth characteristics 
of birds on three commercial farms (Liljedahl, 1970). The 
use of growth equations in describing growth of birds was 


Shown by Ricklefs (1968). 


The basic biological parameters that growth curves can 
deduce are size and rate parameters (Fitzhugh, 1976). He 
stated that the first parameter (A) establishes the position 
of the individual or group in the general size space at a 
given reference age, usually maturity. The second parameter 
(kK) concerns growth rate to body size. When the size 
parameter refers to mature size this parameter defines 
average maturing rate (Fitzhugh, 1976). A third group of 
parameters referred to as the inflection parameters are 
often determined in studies. The point of inflection, or 
the inflection parameter as it is sometimes called, is that 
point at which growth rate is maximum (Brody, 1945; Brown et 
alia, 1976) « Weight at inflection (Y*) and time at 
inflection (t*) have been determined in several studies 
(Timon and FEisen, 1969; Bandy et al., 1970; Brown et al., 


1976) « 


In addition other parameters such as the predicted 
absolute and relative growth rates and weighted absolute 
growth rates can be computed from the functions when the 


values for (A), (B), (k) and (M) are known (Richards, 1959; 
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Fitzhugh, 1976). On this basis in addition to comparing 
fitted growth constants already included in the function, 
derived growth rates can be compared between breed groups or 


populations (Eisen et al., 1969; Rutledge et al., 1972). 


Brown et ale, 1(1972a) stressed the importance of 
providing a true aduit weight in the data as the value of 
(A) is underestimated whenever weights are not available up 
to full maturity. Thus, the fit of the model and derived 
parameters will be determined by the last weight observed in 
the data (Brown, 1970; Fitzhugh, 1976). Furthermore, when a 
growth curve approach is used to study the efficiency of 
selection for growth traits the entire growth curve should 


be considered rather than smali segments (Brown, 1970). 


All growth curves when applied to any form of 
fluctuating data tend to smooth out the irregularities 
(Brown et ale, 1972a; Fitzhugh, 1976; Brown et al-, 1976). 
Brown (1972a) reports that as fluctuations occur when field 
data are used a correction can be made for factors such as 
weight loss during parturition or lactation but difficulty 
arises in adjusting for the corrections when trying to 
interpret the data. Fitzhugh (1976) reported that 
adjustments for weight losses of cows after parturition 
overestimated the mature weight. Fitzhugh (19/76) therefore 
Suggested other alternatives such as skipping the 


overfluctuating data and using weights which conformed 
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closely to the prior weight in relation to a curve. This 
procedure however will lead to some subjectivity in the 
data. Brown (1970) as reported by Fitzhugh (1976) stated 
that simultaneous adjustment of data while fitting non 
linear models was computationally intractable and fitting 


models to unadjusted data seemed most appropriate. 


The reports appearing in the literature utilizing a 
growth curve approach and the types of study reported are 
summarised in Table 1. Most authors cited have worked with 
growth functions that are of particular interest in this 
report, some restricting themselves to single functions 
while Brown et al. (1976) worked with all five functions of 


interest to beef production. 
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studies from the literature 


ee Ee A Species Punction(s) Type of Study 
Brody (1945) Hany Brody General application and 
derivation 
Laird et al (1965) Mice Gompertz Derivation, fitting and 
application 
Laird (1965) Man/mice Gompertz Application of specific growth 
rate to whole and part of peek 
Laird (1966) Birds/mammals Gompertz Application to embryonic 
growth 
Laird (1967) fan Gompertz Application and discussion of 
the human growth curve 
Ricklefs (1967) Birds Brody Application and evaluation of 
‘Logistic grovth rates in migratory 
Gompertz birds 
SSSR CH SHS OS SHO SSS SSCS S5S S56 as eS ASE te en rw ow oo a wo ee ee on sne 
Kidwell et al (1969) Bice Gompertz Analysis of hybridity between 
crosses of inbred lines 
Timon et al. (1969) Hice Richards Fitting, h?, 
Logistic correlations between and 
within parameters and functions 
Eisen et al Mice Bertalanffy Pitting and comparison of 
(1969) Richards the fit of each function 
Logistic 
Gompertz 
Rutledge et al. (1972) Mice Richards Application, fitting ana 
Logistic correlations between parameters 
and functions in a cross 
fostering study 
Eisen (1976) Mice All General review ’ 
Joandet et al Beef cattle Goa pertz Application, fitting, economics 
(1969) 
Fox (1971) Ewes Bertalanfry Application, fitting, h?, and 
Richards correlations between parameters 
Logistic 
Gompertz 
Fitzhugh (1976) Many but with All A review of growth curve analysis 
special reference and an alternative approach to 
to beef cattle study growth 
Brown (1970) Beef cattle Brody. Fitting, 
Bertalanffy application to cattle growth, 
Richards correlation between parameter 
Logistic estimates 
Gompertz 
Brown et al (1972) Beef cattle Brody Application, fitting, 
discussion on maturing rates 
between breeds and sexes 
Brown et al (1976) Beef eattie Brody Application, fitting, discussion 
Bertalanffy of parameter estimates. 
Richards Correlations between paraneters 
Logistic and between functions 
Gonpertz 
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The heritability estimates of the fitted parameters (A) 
and (k) are shown in Table 2, listed by author, species and 


function. 


The heritability of the asymptotic weight appears to be 
high in the studies reported by Timon and Eisen (1969) and 
Fox (1970) with the exception of an estimate based on the 
Richards function in the latter study. Lower heritabilities 
were reported by Brown et al. (1972a). Also Brinks et al. 
(1964) and Calo et al. (1973) have reported higher . 


heritability values for mature weight with actual data. 


The parameter (k) which is the maturing rate showed a 
range in the heritability values from 0.30-0./6. There 
appears to be a difference between the heritability values 
computed from the Richards and Logistic functions (Timon and 
Eisen, 1969). This difference was attributed to the 
increased individual variation in the parameters of the 
Richards function (Timon and Eisen, 1969). In the same 
study genetic and phenotypic correlations between functions, 
in relation to the same parameter were computed and the 
genetic correlation between the (k) values estimated by the 
Richards and Logistic functions was 0.0, while the 


correlation with respect to (A) was 0.98. Timon and Eisen 
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(1969) suggest that the two functions may be estimating 


different traits. 


The inflection point in cattie has been estimated as 
between 6-18 months of age (Fitzhugh, 1976). Timon and 
Fisen (1969) reported heritability values approximating 1.0 
for the age at the point of inflection (t*) and suggest that 
most of the variation is thus of an additive type. Rutledge 
et al. (1972) reported heritabilities of 0.08 and 0.05 for 
(t*) based on the Logistic and Richards functions. The 
authors attributed the low values largely to post-natal 
maternal influences. Heritability estimates of 0.86 for the 
weight at the point of inflection (Y*) have been observed in 


Mice based on the Richards function (Timon and Eisen, 1969). 


The heritability value of the parameter (B), an 
integration constant, is high and Brown et al. (197 2a) 
suggested that early weight changes were thus highly 


heritable. 


Two forms of correlation studies are reported in the 
literature, one restricted to correlations between parameter 
estimates based on the same function (Eisen et ale, 1969), 
and the other involving correlations between functions with 
respect to the same parameter (Brown et al., 1976; Timon and 


Eisen, 1969; Rutledge et al., 1972). 


Correlations between (A) and (k) are shown in Table 3. 
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All correlations, genetic, phenotypic and residual between 
(A) and (k) were negative, ranging from -0.20 to <-0.95 
indicating that animals with higher maturing rates (k) 
reached lower mature weights A) (asymptotes). This 
biological phenomenon is characteristic among species and 
growth functions. Fitzhugh and Taylor (1971) in an 
alternate approach to the study of growth have arrived at 
Similar conclusions about the relationship between maturing 


rates and mature weight in beef cattle. 


In the studies reported by Brown et alj (1972a) and 
Eisen et ale (1969) there appears to be a difference in the 
relationship between (A) and (k) with respect to breeds or 
lines. In the study reported by Eisen et ale (1969) three 
lines of mice were studied and the genetic correlations 
between (A) and (k) were low in the first, medium in the 
second and high in the third line. Brown et al. (1976) 
recognised some degree of independence of the parameters {A) 
and  {k) within breeds and suggests that this genetic 
antagonism can be partially overcome by cross breeding. In 
the report by Brown et al. (1976) six breed groups were 
studied and although the Richards (kK) value was 8/74 greater 
in one breed group compared to the other, there was only a 
2% decrease in the asymptotic weights. Brown et al. (1976) 
suggested that due to the variation in the relationship 


between breeds, the rate of maturing (k) may be amenable to 
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genetic change. 


Correlations between parameters other than (A) with (k) 
from four studies are presented in Table 4, classified by 


species and function. 


Correlations between the inflection parameters (Y¥*) and 
(t*) and size parameter (A) were reported in several studies 
(Timon and EFisen, 1969; Rutledge, 1972; Brown et al., 1976). 
The phenotypic correlations between (A) and (t*) based on 
the Logistic functions (Eisen et al., 1969; Timon and Eisen, 
1969; Rutledge et al., 1972) were positive ranging fron 
approximately 0.3-0.6. Phenotypic correlations between the 
same two traits based on the Richards function were negative 
with values of <-0.10 (Timon and Eisen, 1969) on mice and - 
0.61 (Brown et al., 1976) on beef cattle. However, in the 
first study (Timon and Eisen, 1969) a positive genetic 
correlation of 0.49 was reported. Positive phenotypic 
correlations between (A) and (Y*) were observed by both 
Brown et al. (1976) and Timon and Eisen (1969) suggesting 
that animals heavier at inflection are also heavier at 
maturity. Furthermore, strong positive correlations between 
body weights at different ages are shown by Fitzhugh and 
Taylor (1971), Brinks et al. (1964) and Smith et al. 


(1976). 


Based on the Logistic function, both genetic and 
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phenotypic correlations between the maturing rate parameters 
(kK) and the time of inflection (t*) were negative, the 
values ranging from -0.46 to -0.8 (Timon and Eisen, 1969; 
Fisen et ale, 1969; Rutledge et al., 1972); however, 
correlations (rg and rp) between the same two functions 
based on the Richards were positive, the values ranging from 
0.23-0.77 (Timon and Eisen, 1969). Conversely negative 
correlations between the two traits observed by Brown et al. 
(1976) among beef cattle, suggested that early maturing 
animals reached their inflection age early. Positive 
correlations between (k) and (Y¥*) were observed by Timon and 
Eisen (1969) whereas there was no apparent correlation 
between the same two parameters in the study by Brown et al. 


(1976) . 


Although among mice there appears to be a _ strong 
phenotypic correlation between the time of inflection (t*) 
and weight at inflection (Y*) (rp=0.85) (Timon and Eisen, 
ibe Bebe) = Brown et al. (1976) reported that in beef cattle 
the weight at inflection showed no relationship to the age 


at inflection. 


In summary, correlations between the size parameter (A) 
and rate parameter (k) and the inflection parameters (t*) 
and (Y*) appear to be specific to each function as each 
equation determines a different inflection point. Within 


each function the point of inflection is determined by the 
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function itself rather than the genotype of the animal 


(Fitzhugh, 1976). 


Genetic and phenotypic correlations between (A) and 
absolute growth rate (AGR) in Table 4 were positive, based 
on the Richards function the values ranged from 0.31-0.42 
(Timon and Eisen, 1969). Similar positive correlations were 
observed by Eisen et al. (1969) and Rutledge et al. 
(1972). Negative correlations were reported by Brown et al. 
(1976) between mature weight (A) and absolute growth rate 
(AGR) at the point of maximum growth ie. inflection. 
Positive phenotypic correlations between {k) and  (AGR) 
suggested that animals with lower maturing rates (ie. late 
maturing) showed slower growth rates (Brown et al., 1976; 
Rutledge et al., 1972). Brown et al. (1976) suggested 
that, based on the Richards function a negative association 
between (A) and (AGR) and a positive association between (k) 
and  (AGR) indicated early maturing animals had larger gains 
at inflection but grew to smaller aduit weights. Taran 
studies the correlations between (A) and (M) were negative 
suggesting that heavier animals at inflection reached 


heavier adult weights (Brown et al., 1976). 


The second type of correlation observed in growth curve 
studies are those between functions within a parameter 
(Timon and Eisen, 1969; Brown et al., 1976). The general 


conclusions are that parameter estimates of functions with 
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Similar biological interpretations are positively 
_ correlated, especially the asymptotes but these estimates 
are often smaller than expected (Fithugh, 1976). Fitzhugh 
(1976) further adds that similar parameters from two models 


may not measure the same biological phenomena. 


2.7.2 Alternate growth analysis 


Fitzhugh (1974, 1976) stated that most growth models 
assumed growth is a uniform monotonic increasing function 
with time. Furthermore, he adds that all growth functions 
are insensitive to irregular spacing of size or weight-age 
points and this determines not only their fit but also the 
accuracy; even when animals are weighed at regular 
intervals, the intervals are irregular with respect to 
physiological age and production status. Fitzhugh (1976) 
stated that when working with larger animals under field 
conditions, fitting an appropriate growth equation smooths 
out the irregularities in the actual curve and thereby 
overides the variation due to factors such as climate, level 
of feeding and nutrition, condition or production level, 
thereby possibly obscuring some of the more important 
genetic and environmental phenomena. With some of these 
concepts in mind, Fitzhugh and Taylor (1971) and Fitzhugh 
(1974, 1976) proposed an alternative two component model, 


referred to as an “equation-free analysis of growth curves", 
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for irregular data which preserves the variation in the 
Original observations. Taylor (1965) stated that it has 
long been recognised there is a strong tendency for animais 
with larger mature weights to take a long time to mature. 

Taylor (1968) deduced a relationship equating the time taken 
to achieve a certain proportion of the mature weight to a 
maturing interval, and a degree of maturity. When weight 
(Y) at some point in time is expressed aS some proportion 
(P) of its mature weight, then the proportion Y/A where (A) 
is the mature weight, is referred to as the degree of 
maturity  (U) (Taylor, 1965). According to Fitzhugh (1976) 
the degree of maturity parameter (U) and the parameter  (k) 
from non linear models are alternate ways of characterising 
“animals for being early or late maturing. The parameter (U) 
in addition does not have the disadvantage of (k) as it can 
be applied to small areas of the growth curve, or to 
economically important stages during growth and need not be 
computed from birth to maturity (Fitzhugh, 1976). The 
method of analysis has been used by Smith et al. {1976a, 


1976b) to study growth traits. 
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3. GROWTH FUNCTIONS 


Table 3 lists the Richards, Brody, Bertlanaffy and 
Logistic growth equations and some of their properties 
considered in this report. The four functions were similar 
in their computational formulae and had been used to 
describe beef cattle growth. The parameters or traits and 


the appropriate codes are presented in Table 6. 


3.1 Interpretation of growth curve parameters 


The fitted parameters obtained from a growth curve 
analysis are CAjiee -(B), (kK) fand (M4). Other derived 
parameters such as the weight at point of inflection (Y*), 
age at point of inflection (t*) degree of maturity (u) and 
absolute growth rate  (AGR) are of interest and can be 


interpreted biologically. 


A major parameter each growth function describes is (A) 
the asymptotic weight which is considered to be a measure of 
mature weight; the value of (A) largely depends on the value 
chosen as the adult or mature body weight in the data. The 
value of (A) is computed as time (t) approaches infinity. 
Brown et al. (1976) stated that "the asymptote more nearly 
represents mature weight at a constant condition relative to 
the individuals norm for body composition under a given 


production environment rather than would a single weight", 
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Table 6. Explanation of symbols 
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Effect Code Units 
Asymptote A kg 
Maturing rate k a/day 
Integration Constant B = 
Shape Parameter M = 
Exponent e = 
Time t days 
Absolute Growth rate dy/dt AGR kg/day 
Weight at Inflection Vs kg 
Time of Inflection Boia days 
Weight uf kg 
Weight Predicted MM kg 
Degree of Maturity U a 
Year YR = 
Age of dam AD months 
Bizteh weighteof calf BW kg 
Mean prediction error MEE ri) 
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or a Characterization of the mature body weight of a breed 


group or population. 


The value (k) is the instantaneous relative growth rate 
computed from birth to weaning (Brody, 1945). It was also 
considered as growth rate relative to mature size and was 
referred to as the maturing rate (Taylor, 1965; Brown et 
alse, 1972ab; Brown et al., 1976). In the present analysis 
the maturing rate was expressed as percent per day, as the 
time intervals were presented on a daily basis in the data. 
The values of (k) based on the same function are comparable 
between breeds although comparisons of the parameter between 
functions have been guestioned (Timon and EFisen, 1969). A 
higher value of (k) computed over equivalent weight-age 
points would indicate a rapid maturing type whereas a lower 
value would indicate a slower maturing type. The parameter 
(k) has been used by Brody (1945) to compare growth patterns 
of different species. The use of it in the study of 
maturing patterns among cattle has been demonstrated by 


Brown et al. (1972a,b). 


(B) is an integration constant adjusting for age and 
weight at birth and is used in the computation of other 
parameters such as absolute growth rate but in itself is not 


biologically interpretable. 


Absolute growth rate (AGR) which is the Slope of the 
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growth curve at a point (t) is the first derivative of each 
function. It can also be calculated using the derived 
parameters (A), (B), (kK) and (M). Unlike in a usual 
analysis of daily gains, in a growth curve study absolute 
growth rate can be mathematically calculated at a point or 
Stage in contrast to a period, or interval. Absolute growth 
rate in the present analysis was calculated as gain in 
Kg/day. Derived life time gains such as absolute growth 
rates, maturing rates and relative growth rates could be 


obtained once the values of the fitted parameters are known. 


The weight at point of inflection (Y*) and the age at 
the point of inflection (t*) can also be calculated on the 
Richards, Logistic and Bertalanffy functions once the vaiues 


of (A), (B), (kK) and (M) are known. 


The ratio of (Y¥) to (A) in relation to time measures 
the degree of maturity which can be expressed as a 
percentage. Thus if, 
Y =Aql +t pe- Kt is a generalized growth curve 


then Y /A=(1 + Be- Kt, =Degree of maturity (J) 


From the preceeding discussion it becomes evident that 
a correct value representing the adult weight is exceedingly 
important as the maturing rate parameter (kK) and _ the 


inflection parameters, are computed in relation to (A). 
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4 MATERIALS AND METHODS 


The 6 breed groups, codes and approximate breed 


composition of each breed group are shown in Table 7. 


The Hybrid population was developed from three 
foundation breeds, namely Angus, Galloway and Charolais. 
This Hybrid line was started in 1961 and has been bred 
essentially as a closed line although some introduction of 
outside breeding continued up to 1970. At present the 
approximate breed composition of the Hybrid herd is 1/3 
Angus, 1/3 Charolais, 1/4 Galloway and the remainder divided 
among several other breeds such as Hereford, Holstein, 
Jersey and Brahman introduced through foundation animals 


that were not pure (Berg, 1971, 1975). 


The Hereford line was established in 1960 and has been 
maintained as a purebred line from its initiation with 


limited continued introduction from outside herds. 


The four populations 100% BR, /5% BR, 50% BR 50%LD and 
15% LB were derived from a line referred to as the 
crossbreds, in which the suitability of various combinations 
of Beef type x Beef type or Beef x Dairy types are being 


investigated. These populations will be collectively 


referred to as the XB groups. 
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The data for the Hybrid and Hereford female populations 
were collected over a period of 10 years from 1962-1971 and 
for the four XB groups data was collected over a 6 year 


period from 1965-71 excluding 1969. 


4.2 Breeding and selection program 


Dams chosen to produce the subsequent generations are 
selected according to their reproductive ability, growth 
performance and general soundness. Selection on the basis 
of reproductive performance is strict. Cows are bred in the 
Summer during a two month period, pregnancy tested during 
the winter and those that are not in calf at this time are 
culled. Calving usually occurs during spring the following 
year and is spread over a 2-2 1/2 month period. Calves are 
weaned at approximately 6 months and placed on a 140 day 
feed test over the winter. The feed test period is usually 
from November to mid March during which time all male calves 


are on full feed and females restricted feed. 


Breeding bulls for populations are selected according 
to their pre-weaning and feedlot gain and about 5-7 bulls 
are used within each of the Hybrid and Hereford populations 
each year, one quarter of those selected being from the 
previous year and three quarters being yearling bulis. Ail 


heifers are exposed to the bulis for the first time when 


they are 15-17 months of age. 
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management 


The feeding and management practices of the breeding 
herds conform to that of commercial operations and the type 
of feed used during the experimental period depended on the 
availability of feed and cost. The winter feeding schedule 
is tabulated in Table 8 (Berg, 1971, 1973). During spring, 
summer and fall all dams were allowed access to free 


grazing. 


4.4 Choice of data 


All cows at Kinsella are weighed at intervals of 


approximately 2 months while heifers are weighed monthly. 


In a preliminary study of the data 19 weight age points 
were plotted on a graph to observe the shapes of the growth 
curves and the general variation between weight age points. 
It was noted that the data was of a fluctuating type with 
much of the variation being observed between summer and 
winter weights. Based on visual examination, 9 of the 19 
points were selected for growth curve analysis. The ages 
selected were: birth (0 days), weaning (6 months), yearling, 
15 month weight, 18 month weight, 20 month weight, 32 month 
weight, 44 month weight and 56 month weight (adult). The 


last four weights were winter weights taken during the month 
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Table 8. Winter feeding schedule, Kinsella Alberta, 
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Year 


Age and Type of 
Animal 


Feeding Schedule 


we oe a mr em OO nO 8 Pee Oe we we ce Owe Gere MOTD DH WODMoOwBaA BODE anReenowenrwDeasaeeouescecececcooean 


Prior to 1963 


1963-1965 


19 €4-65 


1965 to 1968 


1968-1969 


1970 on 
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Cows and heifers 


Cows and heifers 


Heifers and 2 year olds 


Heifers and 2 year olds 


Cows and heifers 


Wet Cows 
Keifers 


Cows and heifers 


Heifers 


Heifers 


Winter grazing only 


Hinter grazing 

Barley soybean pellet containing 15% 
CP and 20,000 IU VIT A/LB, fed at 

7 lb/head twice a week 


Additional hay 640 lb/heifer over 
wintering period standard 


Hay ration increased to 1300 lb/head 
Over wintering period due to 
severity of winter 


72 lb bay per week/head on Monday, 
Wednesday and Friday and a barley 
soybean pellet containing 29% CP 
and 70,000 IU Vit A/lb. Additional 
27 ib rolled oats at 9 lb/day 
offered on alternate days Tuesday, 
Thursday and Saturday 


In addition to above 5 lb straw on 
grain days and Sundays or when 
temperature was below 0°F 


75% Straw replaced hay due to cost of hay. 


Supplemental barley soybean fed as usual. 
Grain mixture containing 1 part barley 
and 3 parts rolled oats fed at 30 

lb/head per week. Hay at 4 lb/day 

fed on grain days or when the 
temperature was below 0°PF 


Heifers given 4 lb/hay/day regardless of 
weather conditions in addition to 
straw and grain as above 


we eee meee qeeorae 


39 


ow wes ds iy wee ng « or hy =e shen phe 


a 
seu beh 


(dein ey: nay oe eaten er op ay ate Peg hiatie At: fe 


za Sates 4 she Lek» 
2a fwd) Naat “ei sith 


ara Nene Py hha Ry y} ih 
“praars> es wg 


hao tat “eter Ae ate 18 


aes¥- 


23 Waste ae De Salsa! 
Led 


. Tau BS PONE TY 


sate ayy A fot any 
VA 4seq é a te owl 

ie Pen eae) ar 

Be ae a ed 

rear. 

an 


Pe! wept hse) = 
hee ‘ae bee 


ee ee 


4Q 


of January of each year. 


The number of animals used in each data set along with 
the humber of observations at each age are shown in Table 9 
for the six lines with unadjusted data and the Hybrid and 


Hereford lines with the adjusted data. 


As the animals used each year were born over a 2-2 1/2 
month period this variation existed at each agew Due to an 
overlap of weights at yearling a weight recorded at 11 
months of age and a 15 month weight were regarded as the 
yearling weight in ali populations. Thus, all breed groups 
had more observations recorded at yearling than at the 
preceeding weaning weight. In the adjusted Hybrid and 
Hereford data sets, some animals had an 18 month weight, 
while others had a 20 month weight whereas in the unadjusted 
data only a single animal weight represented 20 months. 
Consequently, more animals were analysed in the adjusted 
data sets. In the XB populations and the adjusted data two 
weights at 18 and 20 months of age were regarded as weight 
representing 20 months. The data analysed was of a 
longitudinal type where a set of observations was available 


on each animal depending on the length of time it remained 


in the herd. 


- a 
, bree ve 
a ue oe 
| rh . 7 
AY EA at, 

(2 Ry 

a we, 7 1 
ee . 
4 , Ui i aM ¥ 7 


iiiy ynoss 264 Atal digg, aiesoetly 


e akvaer ax! eeode ‘oak paal doe pase 


Mei 
eS 5386 odeut be sa a 


mat 4 wh 


ie *, 
i 


th gs Bo pa@oS es yatinesy tg: 


Ck a 


sat Be he Taw, Aus | ekew dsoon: & 5 | 
| ates , set, 


tquosp Beoad ite yk’ : -anaid atagod, Sas, 
wae +e “pad? chara yao pioder _ post 
bie BEV bode rehe! add aes ‘, yu 
yr aso Sree st Ur at bed pene fas 


1 La Poth ce 


oz Hi bom US te on: ants? +" bapaen oat ie 
/ pegeupbs saat ag. apes eee | 
oud ie Be bus aur tig. out Se 


Ae Can, : 


tiview e@F 


# FO BSW 


cay taht aie “ iw: 


ne ja 8 on 


neni asm nou 


41 


DSS SCONe SSS bee SG SeSes eG FS OFS PSSPS S OFELSSSWSASS POOLS MP SHS OSHS SSS SS SSSFOSSSSSOOSSHR SESS OS OSSESSE OS SR SSESSS SSESS 22 See 


6LE ZOL Ent ZBL Ly 6L9 LOE 6LE PIOFOISH 
nSh SSL eee 692 L6S 9LE6 6Eh nS praq ky ‘ peasnfpy 
LZ L OL BL BC ce Lé LZ aT “S2 
8c 9 EL 6L te Oh Be 8c GT %0S Ya %0S 
9S cL a4 6€ OL €6 9S 9S Yd SSL 
He GE Eb Sc Sh cs ne te ud ZOOL 
pth L9 68 £OL 9eb Ly a4 hit fob pirozorsH 
£02 BEL an’ BSL LLL 86E €0z £072 ptaq xy peysnfpeun 


S28 SOB ST8SSSWOS SS SSSS SS &S SOS OS SSSSSBSSVSSSOSST SSOQeIsSe sees seer eseoees QO ee oh es eee’ ewS Sooo ase eso orceooceceoe2 22 weweeos Sa Ge eo oe 


sjTeatue Jo *uow 9¢ “uOw hh  cuom zE *uom OZ  °TIk uvoR q31T@ dnoz6 paeig Peg 


suoT#RAZTeSqO JO TaquNN 


S28 S28 SSS SES ® @ SSS B® FSS OSS VPSSSSB®SSSOeOoee DSW®se®Zeoe 2S ee2 S©SSS € O62 SS SS SSESE SSS OS SSSSSSSTS SESS SS SESS SESS SS SHE Se Tee 


dnoi6 abe pue pasrq yovra ut pasn sTemtue Jo JequnN °6 aTqeye 


42 


4.5.1 Unadjusted 


Growth curves were fitted to a sample set unadjusted 
weight data from of 203 Hybrid and 144 Hereford dams. 
Similarly, growth curves were fitted to unadjusted weight 


age data the four XB groups using all available records. 


4.5.2 Adjusted 


Growth curves were fitted to adjusted data from the 
Hybrid and Hereford lines, the data adjusted for year 
effects and age of dam. All available records comprising 
454 Hybrid and 319 Hereford females were used for the 


analysis. 


Weights of the HY and HE breed groups from birth to 56 
months were adjusted for year (YR).j~ Weights from birth to 
18 months were in addition adjusted for age of dan (AD) 
Classified as 2, 3, 4, 5, 6-8 and >9 years. The model used 
was: 

Wijk = OWL ijky = 701 5 Cy 
where, 

Wijk = adjusted weight 

Wijk = the unadjusted weight of the ith year, 

jth age of dam and kth observation 

Ci = least squares constants for years 


Cj = least squares constants for age of dam 


‘Ipdaaeetionw aba oiijeoe! BUR | 
aii he Soo aig! wef ‘oe! | ie, 3 


Ney ine penne uy | ‘og 7 ‘ 
Adrwondn 1deh tee ee nade 


ee ~~ 


; - 
f x ip. , ; i My iat 
- A a wee A ; My : 
; fa rf Be f me a 
rt ) Set 8 A Nney ' b : o 
op iw Pent My), 
: i \ 


Kad MOAY BAS Ba JaAt hb pee 5a 
eke, | batagy De is Lob . BAe) lat 

1 ' Vi ; ! : ay wt ; 
“‘POLSZLDIBOD ah re $3) ar di ii 5, asi 


ans tok | Bea ce hd 


ae re OF oe, aT 
ae gh ae woz tek 
te mee de pale! i Af 
as fetlon’ gi’? “Et | 


ar bate eed Boa ih | vale ie ie n 
! 7 ee on = son et a 


U 
st : ay 


| . Ww 


dy \ : Pate” 7 a a - ae A 


43 


(Birth to 18 months) 
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The fit of the Richards, Brody, Bertalanffy and 
Logistic equations, over all ages and at 7 weights which 
included birth, weaning, yearling, 20 months, 32 months, 44 
months and 56 months (adult) were studied in detail using 


unadjusted and adjusted data. 


4.6.1 Construction of growth curves 
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The four growth functions - the Richards, Brody, 
Bertalanffy and Logistic were fitted independently to each 
data set. A non-linear least squares programme BMDO7R (BMD 
Health Services Computing Faculty UCLA) which used the Gauss 
Newton iterative procedure fitted the models to the data. 
The programme allowed a maximum of 100 iterations in the 
computation of each set of parameter estimates and was run 


With no constraints. 


4.6.2 Estimation of overall fit 


The accuracy of overall fitting was determined using 
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the square of the multiple correlation coefficient (R2). 


Because in non linear regression the deviation of 
observed minus predicted do not necessarily equal zero and 
because it was expected that at a particular age an equation 
may over or under fit, a measure of average relative 
deviation called mean prediction error (MPE) was developed. 


1@\ 
% MPE = s(yi-?i/yi x 100)/n 
vey 


A lower MPE indicated a lower non symmetrical error or 


a symmetrical error of any size and a better relative 


prediction of weight. 


Overall fit was also compared by calculating MPE over 
all 7 ages and by use of residual variance from a one way 


analysis of variance where differences between ages were 


removed. 
4.6.3 Comparison of fit within age and 


The criterion used to determine the degree of fit of 


each function at each of the 7 ages was the mean prediction 


error estimate. 


The MPE estimates computed from each growth function 
within an age class and breed group were compared by a ‘*t! 


test to establish significant differences (Snedecor, 1962). 
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4.6.4 Comparison of fit between 


unadjusted and adjusted data 


The mean prediction error estimates and residual 
variances were used to compare the fit of each growth 
function between the unadjusted and adjusted data over all 
ages from birth to 56 months and at 7 ages in the Hybrid and 
Hereford lines, The MPE based on each function were 


compared by a 't' test (Snedecor, 1962). 


4.6.5 Analysis of the consistency of fit 
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The overall consistency of fit between breeds was 
determined by comparing the mean prediction error estimates 
within the same function. MPE estimates that were small and 
not significantly different between breeds provided a more 
consistent overall fit. A ‘'t* test (Snedecor, 1962) was 


used to establish significance between means. 


A Duncan's range test was applied to the mean 
prediction error estimates within each function for ail sets 
of data to determine the uniformity of fit. The test was 
made to determine the significant differences among the 


estimated mean prediction errors over the 7 age periods. 
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A random sample consisting of 50 Hybrids and Herefords 
each having a complete set of adjusted weights were used to 
calculate the correlations between (A), (B), (kK), (M), (Y*), 
absolute growth rate at 6 (AGR6), 12 (AGR12) and 18 months 
(AGR18) based on the Richards function and parameters (A), 
(B), (K), (AGR6), (AGR12) and (AGR18) based on the Brody 
function. In calculating the correlations between the 
weight at point of inflection (Y*) and all cther parameters 
both fitted and derived 25 observations per set of data was 


used. 
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The Richards, Brody, Bertalanffy and Logistic functions 
were fitted to the adjusted weight data from the Hybrid and 
Hereford lines. Weights were predicted at birth, weaning, 
yearling, 18 months and 56 months. Absolute growth 
rates/day (average daily gain) were calculated for the 
periods from birth-weaning, weaning-yearling, yearling-18 
months and 18 months-56 months as a difference in weights 
divided by the time interval. 


AGR/day = Y2-¥1/t2-t1 


The observed gain per day and predicted gain/day based 


on each function were next compared within period using a 
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procedure by Gill (1977) designed to make a number of 
independent comparisons (Dunnett type) between means taking 


into account the heterogeneity of the variances. 


SED = / Se/ni + S2/ni 
1 2 
MSD = t 0.01(SED) 


where, 
SED was the standard error of a difference 
MSD was the minimum significant difference 


S2@ and S2 the variance of the observed and predicted 
i 2 


hi the number of observations 
The MSD was compared with the difference between the 
A 
observed and predicted means (Y¥-Y) to establish 


Significance. 


The same statistical procedure was used to compare the 
observed and predicted mean gain over period based on the 
four growth functions for the adjusted data from the Hybrid 


and Hereford breed groups. 
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5. OBJECTIVES OF THE STUDY 
The objectives of the study were to: 


(1) evaluate the accuracy and consistency with which 
the Richards, Brody, Bertalanffy and Logistic functions 
predicted weight (degree of fit) over all ages and at seven 
ages from birth to maturity for several breed groups and 


data sets. 


(2) determine the effect on degree of fit overall ages 
and at each age when data were adjusted for year and age of 


dam effects in the HY and HE breed groups. 


(3) compare asymptotes (A), maturing patterns (k) and 


inflection parameters (Y*) and (t*) among the breed groups. 


(4) evaluate the correlations between fitted and 
derived growth parameters using information from the two 


functions of best fit to estimate responses to selection. 


(5) evaluate how accurately each growth function 


approximates absolute growth rates Or average daily gains. 
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6-2 RESULTS AND DISCUSSION 


6-1 Least squares estimates for adjusted data 


oS a> ataramacs == =, ease Sa 


The least squares constants used for adjusting the 
Hybrid and Hereford weight data for year and age of dam 
effects are shown in Tables 10 and 11. The year effects 
were variable and weights were adjusted accordingly. 
Adjusting for year removed part of the genetic and 
environmental trends. The results showed a gradual increase 


in the age effect as the age of dam increased. 


Cumulative ‘sigmoid* growth curves characteristic of 
each function fitted to the 8 sets of data are shown in 
Figures 1-8. The letters U and A refer to unadjusted and 
adjusted data respectively. The shapes of the Richard's and 
Brody curves were similar in all of the data _ sets. 
Furthermore, the two individual curves could not be 
separately identified in the 100% BR, 50% BR 504LD, 75% LB 


and Hybrid A sets due to considerable overlap. 


The Logistic function overestimated birth weight and 
underestimated adult weight and this feature was 
characteristic across all sets of data. The Bertalanffy and 
Logistic functions also showed a rapid convergence pattern 


to adult weight starting at approximately 750 days of age. 
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The proportions of the total variance accounted for by 
regression in each function and residual variances are shown 
in Table 12 by breed group. The values reported are 
therefore estimates of the accuracy of the iterative fit of 
the parameters and reflect the overall comparative fit of 
the models. In all data sets the Richards and Brody 
functions provided better fits than the Bertalanffy and 
Logistic functions. The Richards, being a four parameter 
model should account for more of the variation than the 
three parameter models. However comparisons of the Richards 
with the Brody model show no improvement in fit in data 


sets. 


6.3.2 Analysis _of the overall fit within breed groups 


The mean prediction estimates, residual mean squares, 
standard deviations, over all ages are shown in Table 13 for 
the unadjusted and adjusted data from the Hybrid and 
Hereford lines and the four XB groups (100% BR, 75% BR, 
50% BR 50% LD and 75% LB). The growth functions applied to 
all of the data sets tended to overestimate weight from 


birth to 56 months of age as all mean prediction error 
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Table 12. Walues of R2 and residual variances by function and breed 


group 
ee n2 eam DP PS PPP BSS MH DH MOST VPHMWSOHMDOMFMO*M@OMBOT PSAYPSTVISG GSH Sem®mevcsseewmececenane 
Breed Group Function & Re Residual Variance 
Hybrid #4 Richards 93.9 4013916.0 
Brody 93.6 4179393.0 
Bertalanffy 92.6 &851124.0 
Logistic 90.6 6186205.0 
Hybrid 2 Richards 93-4 4096481.0 
Brody 93.4 8096726.0 
Bertalanffy 93.0 4355497. 0 
Logistic 9146 5272479.0 
Hereford 3 Richards 94,8 2030715.0 
Brody 94.8 20 45383.0 
Bertalanffy 94.1 2319372.0 
Logistic 92.2 3046321.0 
Hereford 2 Richards S462 2178 768.0 
Brody 941 2209446.0 
Bertalanffy 94.1 2238985. 0 
Logistic 92.8 2706326.0 
100% BR Richards 94.8 499724.0 
Brody 94.8 499927.0 
Bertalanffy 94.4 534475.0 
Logistic 92.9 682313.0 
75% BR Richards 93.6 94657120 
Brody 93.6 950031.0 
Bertalanffy 93.4 984083.0 
Logistic 92.0 11914925.0 
50% BR 50% LD Richards 94.0 509425.0 
Brody 94.0 509454.0 
Bertalanffy 93.6 5 42974.0 
Logistic 92.1 674866.0 
75% LB Richards 94.4 3 38389.0 
Brody 94.4 338826.0 
Bertalanffy 93.8 372223.0 
Logistic 92-0 482848.0 
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Table 13. sean prediction error estimates, residual mean squares, and standard 
deviations by function and breed group 


Standard 
Mean Prediction Residual Deviation 
Breed Group Function error (¢%) Mean Square (%) 
Hybrid 4 Richards -1.287a 145.083 12.363 
Brody -~3.208b 150.950 14.094 
Bertalanffy “5.484 161.832 19.893 
Logistic 17.3334 274.119 47.982 
Hereford & Richards -1,.224a 158.532 12.845 
Brody -1.940a 159.385 12.832 
Bertalanffy -4.125b 165.564 16.965 
Logistic -15.670c 237.675 42.499 
Hybrid 2 Richards -1.287a 144.091 12.341 
Brody -1.21ia 143.958 12.342 
Bertalanffy -3.817b 150.991 156945 
Logistic =14.,320c 205. 126 40.312 
Hereford 2 Richards -1.382b 141,418 12.121 
Brody -0,149a 140.263 413.027 
Bertalanffy ~2.658c 145.077 13.501 
Logistic -12.0454 193.612 34.051 
100% BR Richards -0.762a 106.436 10.655 
Brody -0.618a 106.088 10.727 
Bertalanffy -2.784b 112.596 12.840 
Logistic -13.086c 159. 686 35.343 
75% BR Richards -1.427ab 189.946 13.946 
Brody -0.917a 186.394 14.137 
Bertalanffy -3.209b 204.796 15.765 
Logistic -13.265c 308,801 36.513 
50% BR 50% LD Richards -1.143a 146.621 12.173 
Brody -1.179a 146.732 12.199 
Bertalanffy -3.298b 144.585 14,232 
Logistic -13.706c 165.610 36.178 
75% LB Richards -0.800a 106.057 10.469 
Brody -1.058a 106.713 10.397 
Bertalanffy -3.306a 113.706 13.446 
Logistic -14.112b 163.091 37.906 


tUnadjusted weight data 
2Adjusted weight data. 


Mean prediction error estimates with the different alphabetic letters within a 
breed are significant P<0.05 
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estimates were negative. The results based on the R2 values 
are similar to the results using a mean prediction error 


estimate when the 3 parameter models were compared. 


Based on the magnitude of the residual mean squares the 
Richards growth function provided the best fit to the 
unadjusted weight data from the Hybrid and Hereford lines, 
50% BR 50% LD and 75% LB breed groups; while the Brody 
function offered the best overall fit to the adjusted weight 
data from the Hybrid and Hereford lines, the 100% BR and 75% 
BR breed groups. The Bertalanffy and the Logistic function 
were consistently poor in predicting weights over the entire 


growth curve. 


Brown et al. (1976) used a Similar approach to 
evaluate the fit of five growth models to beef cattle data 
from two locations. According to their results the Brody 
function showed a low residual variance and hence a better 
fit to one set of data while the Richards provided a better 
fit to the other set of data. The Logistic function showed 
the highest residual mean square and hence the poorest fit 
in both locations. The results of the present analysis show 
that the Richards and Brody functions had a better fit than 


the Bertalanffy and Logistic functions. 


Comparing the overall fit between functions within 


breed groups, a lower residual mean square was always 
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associated with a lower mean prediction error. The Richards 
function had the lowest mean prediction error when applied 
to the unadjusted data from the Hybrid and Hereford lines 
and this value was significantly different for the Hybrids 
but not for the Herefords. The Bertalanffy and Logistic 


functions gave poorer fits. 


In the adjusted data from the Hybrid line low mean 
prediction error estimates were obtained for weights based 
on the Brody function followed by the Richards, Bertalanffy 
and Logistic functions. Thus, in the adjusted data the 
Brody function provided the best fit over the entire growth 
curve although, the mean prediction error estimates were not 
Significantly different from the mean prediction error 
estimates based on the Richards function. In the adjusted 
weight data from the Hereford line, the Brody function 
provided the lowest prediction error which was significantly 
different from the prediction error estimate based on the 
Richards function. The Bertalanffy and Logistic functions 
provided a poorer fit and the error estimates were 


significantly higher compared to the Richards and Brody 


functions. 


In the 100% BR breed group the Brody function provided 
a good fit having a small mean prediction error although the 
error estimates based on the Richards was not significantly 


different. Similar trends were established in the other 
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breed groups. In general the fit of the Brody and Richards 
functions were similar and always significantly different 
from the Logistic. The mean prediction error based on the 
Bertalanffy was similar to the Richards in the 75% BR group 
and similar to the Richards and Brody in the 75% LB breed 


group. 


The value of an empirical technique in describing 
growth data depends on the accuracy with which each data 
point is predicted. Brown et al (1976) stated that the 
usual test of the goodness of fit may be inaccurate because 
of the correlated errors among longtitudinal data points. 
Thus, the authors claim that in their study the evaluation 
of fit was subjective, in which the observed and predicted 
weights at each weight-age point were compared. The 
Similarity between the methods using residual mean squares 
and standardised mean prediction errors in evaluating the 
fit of growth functions justifies the latter as a suitable 
procedure. The technique of using prediction error 
estimates in growth curve analysis has not been previously 


reported in the literature. 


6.3.3 Comparison of overall fit _between_ unadjusted 
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estimates based on each growth function for the unadjusted 
and adjusted weight data, and the difference between 
estimates within function and breed group. The overall mean 
prediction error estimates based on the Richards function 
were not significantly different between the unadjusted and 
adjusted data in the Hybrid and Hereford breed groups. The 
Richards function further fitted the unadjusted and adjusted 


data with the same mean prediction error. 


The £1e of the Brody, Bertalanffy and Logistic 
functions was significantly improved by fitting adjusted 


weights in both breed groups. 


The residual mean squares from a one way anova for the 
unadjusted and adjusted weight data are shown in Table 15 by 
function and breed group. Examination of the residual mean 
squares in both breed groups suggests that adjusting for 
year and age of dam effects reduced the residual variance 


within each function. 


In fitting growth models to adjusted weight data one of 
two things could happen. Firstly, the mean prediction error 
estimates could be reduced and secondly, the variation in 
the data sets could be reduced. In all but the Richards 
function fitted to the Hybrid and Hereford data adjusted for 
year and age of dam, the mean prediction error estimates and 


the residual variance were reduced suggesting an improved 
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fit of the models. Although the overall fit of the four 
models was improved by adjusting the data the overall mean 
prediction error estimates based on the Logistic function 
being high, adjusting the data had its greatest practical 
implications in improving the fit of the Brody and possibly 


the Bertalanffy functions. 


6.3.4 Analysis_of the overall consistency of fit 


SS ae 


Table 16 shows the overall mean prediction error 
estimates based on the four growth functions and the 
differences between estimates compared between breeds and 


within each type of data. 


The Richards function provided the most consistent 
overall fit to the two breed groups within each data type as 
none of the breed differences were significant, the mean 
prediction error estimates being small and similar. The 
Brody, Bertalanffy and Logistic functions were less 
consistent as mean prediction error estimates between breed 
groups were different. Although the difference between the 
mean prediction error estimates based on the Bertalanffy and 
Logistic functions fitted to the unadjusted Hybrid and 


Hereford data were consistent between breeds, they were 


large. 
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Table 17. Mean prediction error estimates by function 
in the XB groups 


ee eooee eo eee ow Cow ee@or ees ee DB eODr@e wa@wwrowrea@n@owmtewaewrenaraesaeewoe Gow @ wo we o 


Mean Prediction Error (%) 
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a-None of the mean prediction errors in columns were significantly 
different P<0.05 
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Table 17 shows the overall mean prediction error 
estimated by each function in the XB groups. All mean 
prediction error estimates were similar within function and 
between breeds. However, only the estimates based on the 


Richards and Brody were small as well as consistent. 


The Richards and Brody functions provided the best 
overall fit to all of the data sets. Bertalanffy and 


Logistic functions in general did not provide a good fit. 


Adjusting the data for year and age of dam reduced the 
residual variance in the two breed groups while significant 
reduction of the mean prediction error estimates were 
observed in the estimates based on the Brody, Bertalanffy 


and Logistic functions. 


The Richards function provided the most consistent 
‘estimate of error between breeds and was thus best suited 
for parameter comparisons. 


6.4 Parameter estimates of growth functions 


The values of the fitted parameters (A), (B), (kK) and 
(4) and standard deviations computed by each function in 


each breed group are shown in Table 18. Due to the variable 
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Table 18. 


Breed Group 


Function 


Asyaptote 
(A) kg 


(B) 


Estisates of fitted parameters by fuaction and breed group 


Maturing Rate 


(k) &% 


(5) 
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Hybrid 2 


Hereford # 


Hereford 2 


100% BR 


75% BR 


50% BR 50% LD 


75% LB 


Richards 
Brody 
Bertalanffy 
Logistic 


Richards 
Brody 
Bertalanffy 
Logistic 


Richards 
Brody 
Bertalanffy 
Logistic 


Richards 
Brody 
Bertalanffy 
Logistic 


Richards 
Brody 
Bertalanffy 
Logistic 


Richards 
Brody 
Bertalanffy 
Logistic 
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Brody 
Bertalanffy 
Logistic 


Richards 
Brody 
Bertalanffy 
Logistic 


tUnadjusted weight data 
Z2Adjusted weight data 


414.65 
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465.09 
$41.07 
425.81 
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458.91 
432.40 
415.8] 
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495.60 
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humber of animals in the sets of data, the parameter 
estimates of the Hybrid and Hereford breed groups will be 
discussed separately from the parameter estimates of the 


four XB groups. 


6.4.1 Asymptotes (A) 


The mean observed and predicted adult weight (56 
months) and the asymptotes estimated by each growth function 


for each data set and breed group are shown in Table 19. 


In the Hybrid and Hereford populations with the 
unadjusted weight data, the Richards function estimated the 
highest asymptote followed by the Brody, Bertalanffy and 
Logistic functions. However, in the adjusted data the 
asymptotes fitted by the Brody function were slightly higher 
than the Richards. The asymptotes fitted by the Brody, 
Bertalanffy and Logistic functions were similar in both data 
sets within breed and function while the asymptotes fitted 
by the Richards function were higher in the unadjusted data 
sets of both breed groups. Furthermore, in the adjusted 
weight data, the asymptotes fitted by the Richards function 
were lower than the observed adult weights, while the 


reverse was true in the unadjusted data sets. 


In the Hybrid breed group using unadjusted data adult 


weight was underestimated by 9.16 kg and 19.72 kg based on 
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the Richards and Brody functions respectively while adult 
weight was underestimated by 19.12 kg and 18.66 kg based on 
the same functions using the adjusted data. Similarly among 
the Herefords the differences between the observed and 
predicted adult weight based on the Richards was greater in 
the adjusted data which reflected the lower fitted 


asymptotes in these data. 


Although all weights were adjusted, adjusting the 
weights at 32 months, 44 months and 56 months of age for the 
year effects in the Hybrid data was more critical in 
determining the asymptote and decreased the mean weight at 
each of these ages (Table 10). Similarly, in the Hereford 
breed group mean 32 month and 44 month weight was decreased 
and 56 month weight increased (Table 11). The general 
reduction of mean body weights during later ages by 
adjusting the data would have accounted for the lower 
asymptotes and mean predicted adult weights in the Hybrid 
breed group. The effects of adjusting were however less 
pronounced in the Hereford breed group. In general 
predicted adult weight was underestimated to a greater 
degree by adjusting the data, which reduced the asymptotic 
weights fitted by the Richards model in particular. The 
asymptotes based on the Brody, Bertalanffy and Logistic 


models were relatively insensitive to the changes in weight 


due to adjusting. 
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Brown et al. (1976) in a study comparing the relative 
efficiency of 5 non linear models in describing growth of 
Herefords fitted asymptotes of 505 kg, 508 kg, 488 kg and 
481 kg using the Richards, Brody, Bertalanffy and Logistic! 
functions respectively, compared with fitted asymptotes of 
476 kg, 464 kg, 433 kg and 420 kg in the present analysis 
using unadjusted weights based on the Richards, Brody, 
Bertalanffy and Logistic functions. The 30-60 kg difference 
between asymptotes of the two studies could be due to the 
choice of weights and the time at which weights were 
recorded in relation to season or to time differences in 


environment. 


In every set of data under study the Bertalanffy and 
Logistic functions fitted asymptotes that were below the 
observed adult weight and almost equal to the predicted 
adult weight. A rapid convergence pattern and 
underestimation of the asymptote is therefore characteristic 
of the two functions. When the estimated shape parameter 
(M4) of the Richards function was close to 1 (Table 18) the 
fit and consequently the asymptotes computed by the Richards 
were almost identical to the Brody, as the computational 


formulae of the two functions were identical when M=1. 
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Table 20 shows the mean observed and predicted adult 
weights and fitted asymptotes by function for the four XB 
groups. In the four XB groups with the exception of the 75% 
LB breed group the Brody function fitted the highest 
asymptote. The Richards function estimated a high asymptote 
to the 75% LB group. However, the asymptotes fitted by the 
Richards and Brody functions were similar between groups and 
differences never exceeded 7 kg. Furthermore, in all XB 
groups the fitted asymptotes were lower than the mean 
observed adult weight and in every instance adult weight was 
underestimated. The Logistic was similiar to the Bertalanffy 
in fitting low asymptotes which closely approximated the 
mean predicted adult weight. Due to the fewer number of 
observations available at 56 months of age and the variation 
in the available observations the fitted asymptotes and 
predicted adult weights may not adequately characterize each 


breed group. 


6.4.2 Maturing rates_(k) 


The parameter (k) is a general maturing rate parameter 


calculated over the entire growth curve period. 


Comparing the (k) values from the Richards function 
(Table 18) using the unadjusted weight data there was a 28% 


difference in maturing rates, Herefords being earlier 
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Table 20. Mean observed, predicted adult weights and asymptotes 
by function for the XB groups 
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Mean Mean 
Observed Predicted 
Breed Group Function Adult wt {kg) Adult wt (kg) A 
56 month 56 month (kg) 
100% BR Richards 482.92 458.49 463.18 
Brody 459.75 465.09 
Bertalanffy 440.46 441.07 
Logistic 425.83 425.81 
75% BR Richards 486.39 448.05 452.01 
Brody 452.76 458.91 
Bertalanffy 431.71 432.40 
Logistic 415.85 415.87 
50% BR 50% LD Richards 530.81 487.94 494.69 
Brody 488.50 495.60 
Bertalanffy 466.47 467.31 
Logistic 452.46 452.47 
15% LB Richards 495.71 473.16 482.79 
Brody 470.88 478.66 
Bertalanffy 450.54 ; 451.56 


Logistic 437.68 437.71 
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maturing than Hybrids, while a 17% difference was observed 
with the adjusted data. On the other hand based on the (kK) 
values derived from the Brody function in the unadjusted 
weight data the Herefords showed a 7% lower maturing rate. 

Brown et al (1976) reported similar results while comparing 
the maturing rates derived by the Brody and Richards 
functions among Jerseys and Herefords. Based on the 
Richards function the Herefords matured later than the 
Jersey while based on the Brody function the Jerseys matured 
jater. However no attempt was made to explain this 


difference. 


Comparisons of maturing rates based on the growth 
models were reported between breed groups (Brown et al, 
1976) and between sexes within breeds (Brown et al, 1972a). 
The values of the maturing rate parameter (k) computed on 
Eten model in the present study were comparable with the 


results in the literature. 


The values of (k) based on the Richards and Logistic 
functions differed widely in each set of data and similar 
observations were reported by Brown et al (1976) and Timon 
and Eisen (1969). In the study of Timon and Eisen (1969) 
heritability estimates for the parameter (k) based on the 
Richards and Logistic functions were 0.30 and 0.76 
respectively. The genetic correlation between the parameter 


(k) calculated from the two functions was 0.0. The authors 
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therefore suggested that the two functions may be estimating 
different parameters, In contrast, the correlations between 
functions of the asymptotes (A) were highly positive and 


approached unity (Timon and Fisen, 1969). 


The differences bet ween the (k) values of the 
unadjusted and adjusted data based on the Richards was in 
part due to the lower asymptotes fitted when the data were 
adjusted. The differences between unadjusted and adjusted 
(k) values based on the Brody function were less pronounced 
between data sets within a breed group because the 
asymptotes were reduced to a lesser degree by adjusting the 


data. 


The Richards and Brody functions were consistent in 
fitting the XB groups with small and uniform mean prediction 
error estimates and were thus ideal for parameter 
comparisons (Table 17). However the Richards function 
provided a poorer fit to the 75% BR group. Comparing the 
values of (k), among the XB groups (Table 18) the lowest 
values of (k) were observed in the 75% LB breed group, 
suggesting that the 75% LB dams matured slower than the 
other three breed groups. British breeds of cattle, namely 
Angus and Hereford are early maturing (Berg and Butterfield, 
1976) and in this analysis an increase in the proportion of 
large beef blood would be expected to reduce maturing rates. 


Based on the Brody function the values of (k) were lower for 
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the 75% BR group ccmpared to the 100% BR group. It is 
therefore indicative that the 75% BR dams matured later than 
the 100% BR dams. The maturing rates computed by all growth 
functions for the 50% BR 50% LD breed group were lower than 
the maturing rates of the 75% BR group suggesting that the 


British/Dairy crosses (50% BR 50% LD) are later maturing. 


In general, the 75% LB breed group showed a late 
maturing pattern followed by the 50% BR 50% LD, 75% BR and 
100% BR groups. The maturing rates in the four XB groups 
appeared to be directly proportional to the amount of 
British beef blood. However, because of the small size of 
Sample analysed the resuits cannot be considered as being 


absolutely conclusive. 


6.4.3 Inflection parameters (Y¥*),_(t*) 


The weight at point of inflection (Y*) and the age at 
which inflection occurs are shown in Table 21 by function 
and breed group for the unadjusted and adjusted data. As 
the Logistic and Bertalanffy functions have inflection 
weights that are fixed (Table 5), there is very little 
variation in the estimates of (Y*) and (t*). The Richards 
which has a variable inflection point allows variation in 
the estimates. Based on the Richards function the mean 


weights at point of inflection (unadjusted and adjusted) 
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were 171 kg and 182 kg for the Hybrid and Hereford groups. 

Fitzhugh (1976) states that as the growth curve is 
essentially linear at inflection the estimated point of 
inflection is influenced more by the type of function used 
rather than the animals genotype. In cattle the time of 
inflection occurs at a point between 6 to 18 months of age 
(Fitzhugh, 1976). The time of inflection for the Richards 


function could not be determined in this study. 


The Brody and Richards function fitted larger 
asymptotes in all breed groups compared to the Bertalanffy 


and Logistic functions. 


All growth functions fitted higher asymptotes to _ the 
Hybrid data compared to the Herefords illustrating that the 


Hybrids achieve higher mature weights. 


In the adjusted data due to a general reduction of body 
weight during the later ages, the asymptotes fitted by the 
Richards function was markedly lower in both breed groups 
while the Brody, Bertalanffy and Logistic functions were 
relatively insensitive to such changes and fitted similar 


asymptotes to the unadjusted and adjusted Hybrid and 


Hereford data. 
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The Bertalanffy and Logistic functions converged 
Sharply at the asymptote therefore resulting in lower 
asymptotes. The predicted adult weight and the asymptotes 
were markedly similar within function. The asymptotes 
fitted by the four growth functions for the Herefords were 


comparable with reports in the literature. 


Based on the Richards function which offered a 
consistent fit to all of the data the Herefords were 
recognised as being early maturing compared to the Hybrids. 
The (k) values fitted by the Richards and Logistic functions 
within breed and data set were different suggesting that the 
two models may be fitting parameters of different biological 


Significance. 


The differences in the (k) values fitted by the 
Richards and Brody functions between the unadjusted and 
adjusted data within breeds are attributed to the 


differences in the fitted asymptotes. 


Based on a limited number of observations the 75% LB 
breed group was later maturing compared to the other XB 
groups. The maturing rates were proportional to the amount 


of British beef blood present in the four XB groups. 
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Table 22 shows the mean prediction error estimates, 
Standard errors, and differences between unadjusted and 
adjusted mean prediction error estimates by function at each 
age for the Hybrid breed group. The minima and maxima for 
the &% mean prediction error by function and age for the 


Hybrid data sets are shown in Appendix Table 1. 


Birth: In the unadjusted Hybrid weight data the 
Richards function showed the smallest mean prediction error 
0.67% for birth weight followed by the Brody, Bertalanffy 
and Logistic functions. The 3 latter functions 
overestimated birth weight; the Logistic function 
overestimating weight by 126.6% The mean prediction error 
estimates computed by the Richards function was 
Significantly less than from the other three. The 
Bertalanffy and Logistic functions predicted minima and 
maxima which were both negative and thus reflected their 
poor performance in predicting birth weight (Appendix Table 


1) 


In the adjusted Hybrid data the Brody function provided 
the best fit to the weight at birth followed by the Richards 


nodel. The mean prediction errors calculated on the Brody 
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and Richards functions were not Significantly different 
although they were significantly smaller than those from the 


Bertalanffy and Logistic functions. 


In general, the Richards function was consistently good 
in describing birth weight in the two sets of Hybrid data 
and the Brody was equally good for the adjusted data. The 
Bertalanffy and Logistic functions were unsatisfactory in 
predicting birth weight as they were associated with large 


mean prediction errors. 


Weaning: At weaning the Brody function provided a good 
fit in both sets of data. The Brody function provided an 
error that was significantly less than the Richards in the 
unadjusted Hybrid data but not significantly different ina 
the adjusted data. In the unadjusted Hybrid data the 
Richards function tended to overestimate weaning weight 
while the Brody, Bertalanffy and Logistic functions 
underestimated weaning weight. However, weaning weight was 
underestimated by only the Bertalanffy and Logistic 
functions in the adjusted data. The Bertalanffy and 
Logistic growth functions provided poor estimates of weaning 
weight as the mean prediction errors were significantly 
larger than those computed on the Richards and Brody 


functions in both sets of data. 


Yearling: At the yearling weight, all mean prediction 
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error estimates computed by each growth function applied to 
the unadjusted Hybrid data set were small. A comparison of 
mean predicticnh error estimates applied to the same 
population showed the mean prediction error estimates for 
the Richards, Brody and Bertalanffy functions were similar, 
and these were significantly different from the error 
estimated by the Logistic function. In the adjusted Hybrid 
data, the magnitude of error was similar although the mean 
prediction error computed by the Bertalanffy function in the 
adjusted Hybrid data was significantly higher than the 
Logistic. Thus, considering both the unadjusted and 
adjusted data the four growth functions gave a satisfactory 
prediction of yearling weight; the Brody function ranking 
the best among the unadjusted yearling weights and the 
Logistic ranking the best among the adjusted yearling 


weights. 


20 months: The weight at 20 months was best predicted 
by the Brody function with less than -0.5% error in the 
unadjusted data and by the Bertalanffy with less than a 
0.61% mean prediction error in the adjusted data. The mean 
prediction error estimated by the Brody function was 
Significantly lower than the mean prediction errors 
estimated by the other growth functions in the unadjusted 
Hybrid data while the mean prediction errors based on the 


Bertalanffy function was significantly lower than the mean 
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prediction errors based on the other growth functions in the 
adjusted data. The mean prediction errors based on the 
Logistic function were consistently larger in all of the 


data sets. 


32___months: At se months there was a general 


overestimation of weight by all four growth functions in 
both sets of data. The mean prediction error estimates were 
also higher in all functions compared to the mean prediction 
error estimates at yearling and 20 months. In the 
unadjusted data the Richards function provided a better fit 
although its mean prediction error was not significantly 
different from that of the Brody function. Similarly, no 
Significance was established between the mean prediction 
errors estimated by the Brody and Bertalanffy functions and 
the Bertalanffy and Logistic functions. In the adjusted 
Hybrid data the Brody function provided a relatively low 
prediction error although the error estimates computed by 
the other functions were not significantly different from 


the Brody function. 


The standard errors of the mean prediction error 
estimates at 32 months were higher in all sets of Hybrid 
data than standard error estimates during all intermediate 
ages. During this time, cows usually lost weight due to an 
imposed temporary winter stress and were ina state of 


negative growth. Furthermore, as animals had passed the 
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point of inflection their growth rates were decreasing. A 
combination of these factors can cause Significant 
deviations from a normal growth pattern and this effect is 
recognised when growth equations are fitted to these types 
of data. AS each growth curve assumes a characteristic 
growth pattern it is not sensitive in reacting to sudden 
fluctuations resulting in a decrease in weight and thus will 
overestimate weight at such stages, as illustrated, at 32 
months of age. In generai, all functions gave poor 
estimations of weight at this age with mean prediction 


errors ranging from 7% to 11.6%. 


44 months: At 44 months of age the Bertalanffy function 
provided the best fit to ali of the Hybrid data. The mean 
prediction error estimates computed by the Bertalanffy 
function were not significantly different from the estimates 
of the Richards, Brody and Logistic functions in the 
unadjusted data but were significantly different from the 
estimates of the Richards and Brody functions in the 
adjusted Hybrid data. All mean prediction errors were small 
and had low standard errors. In general, the weight at 44 
months was predicted well by all growth functions with the 


Bertalanffy function having the smallest mean prediction 


errors. 


56 months: At 56 months (i¢. adult) the Richards 


function provided the best prediction of weight in the 
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unadjusted data followed by the Brody, Bertalanffy and 
Logistic. The mean prediction error estimated by the 
Richards function was significantly lower than those from 
the Bertalanffy and Logistic functions. In the adjusted 
Hybrid weight data set, the Brody function ranked first 
followed by the Richards, Bertalanffy and Logistic 
functions. However the difference between the Brody and 
Richards were not significant. The Bertalanffy and Logistic 
functions provided a poorer fit in both populations because 


adult weight was underestimated (Table 19). 


6.5.2 Comparison of fit between unadjusted and adjusted 


<> ee ee > eS ee ee eS Se Se =n @E a a eee —o_ wea as 


weight data_ within function and age 


The differences between the unadjusted and adjusted 
mean prediction error estimates by age and function for the 


Hybrids are shown in Table 22. 


The birth weights predicted by the Richards function 
for the unadjusted and adjusted data were similar in that 
both mean prediction error estimates were of the same 
Magnitude but different only in sign. Adjusting the data 
significantly improved the fit of the Brody function and in 
fact provided a lower mean prediction error estimate than 
the Richards. The mean prediction error estimates computed 


by the Bertalanffy and Logistic growth functions were also 
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improved by adjusting the data. 


The prediction of weaning weight based on the Richards 
and Brody functions were significantly improved by adjusting 
the data. No significant reduction in the mean prediction 
error estimates were observed by adjusting weaning weights, 


based on the Bertalanftfy and Logistic functions. 


The prediction of yearling weight was not improved by 
adjusting the data. The mean prediction error estimates for 
the unadjusted data based on the four growth functions were 
Significantly lower than the mean prediction errors based on 


the same growth functions with the adjusted data. 


At 20 months of age no significant differences were 
observed between the mean prediction error estimates based 
on the Richards function between the unadjusted and adjusted 
weights. The mean prediction error based on the Richards 
applied to the adjusted data was very Similar to the mean 
prediction error computed on the same function in the 
unadjusted data. The mean prediction error estimated on the 
Brody function was significantly lower in the unadjusted 
data. The Bertalanffy and Logistic functions however 
provided a better fit to the adjusted data at 20 months of 


age and the differences were significant. 


At 32 months of age, no significant differences between 


mean prediction errors of the unadjusted and adjusted data 
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were observed. All error estimates were large during this 
age. The mean prediction errors estimated for the adjusted 
data based on the Logistic function were lower when compared 
to the mean prediction error estimate in the unadjusted 


weight data, but the differences were not significant. 


Mean prediction error estimates based on the four 
functions and applied to the unadjusted and adjusted weight 


data were similar at 44 months of age. 


At 56 months of age mean prediction errors based on the 
Richards function were significantly higher for the adjusted 
data. However no significance was observed between mean 
prediction error estimates based on the Brody, Bertalanffy 


and Logistic functions between unadjusted and adjusted data. 


6.523 Comparison of the variation between unadjusted 
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The standard deviations of the mean prediction error 
estimates within function and breed group for the unadjusted 


and adjusted data are given in Tabie 23. 


Adjusting the data for year and age of dam reduced the 
variation in the mean prediction error estimate for birth 


weight in all growth functions thereby improving the fit. 


At weaning, variation in the mean prediction error was 
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reduced by adjusting the data in all but the Bertalanffy 
function but the improved fit was only observed with the 
Richards. No reduction in the variance by adjusting was 


observed at yearling age. 


The standard deviation of the mean prediction error 
estimates of ali functions at 20 months of age were lower 
for the unadjusted data compared with the adjusted data. 
Hence, the predicted 20 month weights Were closer to the 


observed in the unadjusted data. 


A reduction in the standard deviation by adjusting the 
data was observed at 32 months for all functions. However, 
weight was poorly described at this age aithough the 
variation in the mean prediction error estimates were 


reduced by adjusting the data. 


The standard deviations of the mean prediction error 
remained almost constant at 44 and 56 months of age in the 


unadjusted and adjusted weight data based on the u 


functions. 


Table 24 shows the mean prediction errors at each age 
derived from the Richards, Brody, Bertalanffy and Logistic 


growth functions for the unadjusted and adjusted data from 
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the Hybrid breed group. 


in the unadjusted data the Richards function showed the 
least significant differences between mean prediction error 
estimates ranging from -6.9 to 2.0 and falling into three 
Significance classes with much overlap. The mean prediction 
errors estimated from the Brody function were spread over a 
wider range and showed greater significant differences 
between ages; means were grouped into four significance 
classes. The mean prediction error estimates of the 
Bertalanffy and Logistic functions fell into six 
significance classes, suggesting that the functions provided 


a less consistent fit from birth to 56 months. 


In the adjusted data the mean prediction errors based 
on the Richards function ranged from -8.4 at 32 months of 
age to 3.0 at 56 months of age falling into four 
significance classes. The mean prediction error estimates 
computed by the Brody function in the adjusted data fell 
into four significance classes, and therefore the Brody was 
as consistent as the Richards function. The Bertalanffy 
function was very inconsistent in predicting small and 
similar mean prediction error estimates; the mean prediction 
error estimates fell into six significance classes. The 
Logistic function followed the same trend as the Bertalanffy 
and in addition, the mean prediction error estimates applied 


to the adjusted data were in six significance classes. 
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The Richards was the only function that provided a good 
and uniform fit to all of the birth weight data among 
Hybrids; the function predicted birth weight with less than 
+ 0.7% mean prediction error in all of the data sets. The 
Brody provided a differential fit at birth providing mean 
prediction errors that were small and similar to the 
Richards in the adjusted data but significantly different 
from the Richards in the unadjusted Hybrid data. The 
Bertalanffy and Logistic functions were very poor in 


describing birth weight in both of the Hybrid data sets. 


Weaning weight was best described by the Brody 
function. Significant differences in the mean prediction 
error estimates were observed between the means computed on 
the Richards and Brody function in the unadjusted data while 
the mean prediction errors were small and similar in the 


adjusted data. Thus, the Richards and Brody functions 


described weaning weight adequately. 


Yearling weight was uniformly described by all four 
growth functions with mean prediction errors that were 


Similar in each set of data. 


The Bertalanffy function provided the lowest estimate 
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of error at 20 months of age in the adjusted data while the 
Brody function provided a good description of weight at 20 
months in the unadjusted data. The Richards, Brody and 
Bertalanffy functions all provided an adequate fit at 20 


months of age. 


Weight at 32 months was overestimated by all growth 
functions as they were insensitive to temporary and sudden 
changes in weight. High standard errors were also 
associated with the mean prediction error estimates during 


this period, 


All growth functions described BG month weight 
adequately. In general, during the intermediate ages 


(yearling to 44 months) all functions fitted the data well. 


Adult weight was best described by the Brody and 


Richards functions poorly described by the Bertalanffy and 


Logistic functions. 


Adjusting the data for year and age of dam 
significantly improved the fit of the Brody, Bertalanffy and 
Logistic functions at birth by reducing the mean prediction 
errors and the variance. The prediction of weight and fit 
of the Richards and Brody functions were significantly 
improved at weaning, by adjusting the Hybrid weight data. 
However, the fit of the functions at the yearling stage was 


significantly better with the unadjusted data. Adjusting 
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the weights after 20 months of age did not significantly 
reduce the mean prediction error estimates with the 


exception of the Logistic function at 32 months. 


Thus, the value of adjusting the Hybrid data was 
limited to improving the fit of the functions during the 
pre-yearling ages. A better fit was often associated with a 
reduction in the variance of the mean prediction error 


estimates. 


The Richards function provided a very uniform 
estimation of mean prediction errors in each Hybrid data 
set. The prediction error estimates ranged within narrow 
limits, were smali in magnitude and fell into a lower number 
of significance classes compared with the other three growth 
functions, The mean prediction errors based on the Brody 
function ranked second with respect to the consistency of 
fit between ages, and the uniformity was pronounced in the 
adjusted data. The Logistic and Bertaianffy functions were 
very inconsistent with variable mean prediction errors at 


each age. 
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6.6 Analysis of Hereford data 


function applied to each age 


Table 25 shows the mean prediction error estimates, 
standard errors and differences between unadjusted and 
adjusted mean prediction error estimates by function at each 


age for the Hereford breed group. 


The minima and maxima for the % mean prediction error 
by function and age for the Hereford data sets are shown in 


Appendix Table 2. 


Birth: The Richards function provided the best 
prediction of birth weight in the unadjusted Hereford data 
much like the observations in the Hybrid breed group (see 
Table 22) and underestimated birth weight slightly while the 
other functions overestimated it. The Bertalanffy and the 
Logistic functions were particularly poor in predicting 
birth weight. The mean prediction error estimates 
calculated on all functions were significantly different 
from each other, In the adjusted Hereford weight data a 
similar trend was recognised with the Richards function 
providing a good prediction of birth weight followed by the 


Brody, Bertalanffy and Logistic functions. 


Weaning: The lowest mean prediction error for Herefords 
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at weaning was obtained for the weights based on the Brody 
function followed by the Richards, Bertalanffy and Logistic 
in the unadjusted data. The mean prediction error estimates 
based on the Brody and Richards functions were similar 
whereas, the Logistic computed an error that was 
Significantly higher than the estimates of the Brody and 
Richards functions. In the adjusted Hereford weight data 
the mean prediction error for the Richards function was 
Significantly lower compared to the other growth functions, 
and was therefore different from the trend in the Hybrid 


breed group (Table 22). 


Yearling: At one year of age all functions described 
growth adequately with a minimum of error in the unadjusted 
data much like the results in the Hybrid breed group. 
Furthermore, as in the Hybrids, the Richards, Brody and 
Bertalanffy functions estimated yearling weight with <1% 
error and were significantly lower than the errors computed 
by the Logistic function. In the adjusted Hereford data all 
growth functions overestimated yearling weight with <2% 


error with the Logistic function providing the lowest mean 


prediction error. 


20 months: At 20 months all functions overestimated 
weight with the Richards providing the lowest mean 
prediction error and the Logistic» the highest: in the 


unadjusted data. Significant differences were observed 
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between the mean prediction errors estimated by the 
Richards, Bertalanffy and Logistic functions while no 
Significance was observed between the mean prediction errors 
of the Richards and Brody functions. Similarly, in the 
adjusted Hereford data the the Richards function provided 
the smallest mean prediction error but the estimate was not 
Significantly different from the means based on the Brody 
and Bertalanffy functions, Thus, the weight at 20 months 
was best described by the Richards function followed by the 


Brody and Bertalanffy functions. 


32 months: Large mean prediction errors and standard 
deviations were observed at 32 months of age in the 
unadjusted and adjusted Hereford data for all functions much 
like the trends observed among the Hybrids, for the same 


reasons as previously discussed. 


44 months: The weight at 44 months of age was best 
described by the Brody and Richards functions followed by 
the Bertalanffy and Logistic functions in the unadjusted 
data. No significance was established between the mean 
prediction error estimates based on the Richards, Brody and 
Bertalanffy functions although the mean prediction error 
estimate based on the Logistic function was significantly 
different from the Richards and Brody functions. In the 


adjusted Hereford data, the fit of the functions was similar 


to the Hybrids (Table 22). 
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26 months: For the most advanced age group (56 months) 
the mean prediction error estimates based on the Richards 
and Brody functions were similar and ranged from 1-2% in the 
unadjusted weight data while the Bertalanffy and Logistic 
functions predicted 56 month weight with a significantly 
higher mean prediction error. Adult weight was 
underestimated by ail functions, In the adjusted data no 
Significance was established between the errors computed by 
the Richards and Brody functions and Richards and 
Bertalanffy functions although these values were 
Significantly different from the mean prediction error 
estimates based on the Logistic function. In general, the 
Richards and Brody function provided a satisfactory 
description of 56 month weight in the unadjusted and 


adjusted data. 


weight data within function and age 
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The differences between the unadjusted and adjusted 
mean prediction error estimates by age and function for the 


Herefords are shown in Table 25. 


The mean prediction error estimates at birth calculated 


on the Brody, Bertalanffy and Logistic functions differed 


significantly between the unadjusted and adjusted Hereford 


Le oy a, 
an hee 
{ , f “ eh 
i { ; ; a aa 
iy _ d 
t ‘ i ti 
mie , 
Wil 
‘ ; i 
‘ Wisi 
’ Ars AL 
Z ' iy io 
j ’ er a 
i f ra ' ' ns ee 
be DS rec ae} ECA: SS. Da * 
q tj ‘ % te 
: (ay 
, a a J 
f R 
mvt bl «epg i ve 
ee Dae Aoxge. a 
a7 he 
Fi Ets TCO anki anh? 
- ¢ ale 
: iy eh peas | 
: ~ il Rahints ee i 
i ” . 
} ; aut } ou ri 
i t 
4 ) ey Sn a 
i Ure eA ite ‘ 
pote 0s, . 
Ae alee dem be ie Ee & ’ 
+ 4 M ’ 
4 i) Lee ‘at a 
i, 7 mes, hs 
; : ay ‘ 
‘ , ¢ q 
OF S28 gr 56 14a Si a pare 
Tes 5 : oa Shot 
pa fe hs ye Ker: 
a 7 


y : ia 
j “ \ : Yeah ' ie yo 
o : a OL, 4's , " 
b ul he ete) ills 
ing tae Ay an Lo 
i } Wy ‘hy i , 
: ’ ’ Los ine x J q } vet rh, H 5 
“ ers an A 
; ha ig f F Pid i é i ¢ ” few . a ‘is ¥ 
_— = ae a ea r ny : yi bo So a ; ti wor, at ae 
Fame, D2), 2 eas eee BE then fore ee ~ ‘ . ed pea a sialic * van ibe 
t 


f iss, i i 
= ny 7 
Rs 2 , 
Lie > a fi 
‘ : om " _ iu f ‘ i ; “a a 


cane a a 
aA ae 


al hee 7 : 


107 


weight data while no significance was established on the 
mean prediction errors based on the Richards function. 

Adjusting the data did not improve the fit of the Brody 
function at birth although the fit of the Bertalanffy and 


Logistic functions was improved. 


At weaning the mean prediction errors based on the 
Richards function were reduced significantly by adjusting 
the data while the reverse occurred with respect to the 


Brody function. 


The differences between mean prediction error estimates 
calculated on the Richards, Brody and Bertalanffy functions 
for the unadjusted and adjusted data, at yearling were not 
Significant. The mean prediction error estimates based on 
the Logistic function were significantly lower in the 
adjusted data compared with the unadjusted data. The growth 
functions in general provided a better fit to the unadjusted 


yearling weights. 


At 20 months of age the mean prediction error estimates 
based on the Richards function between the unadjusted and 
adjusted data were not significantly different although the 
Brody function provided a better fit to the unadjusted data. 
The fit of the Bertalanffy and Logistic functions was 


significantly improved by adjusting the data. 


For all other ages, the mean prediction error estimates 
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of the unadjusted and adjusted data within function were 
Similar with the exception of the mean prediction errors 
based on the Richards function at 56 months of age where a 


better fit to the unadjusted Hereford data was observed. 


6.6.3 Comparison of the variation between unadjusted 
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Table 26 shows the standard deviations of mean 
prediction errors for the unadjusted and adjusted data by 
function and age for the Hereford breed group. The Brody, 
Bertalanffy and Logistic functions fitted to the adjusted 
data reduced the variation in the mean prediction error 
estimates at birth and a reduced variation was associated 
with a closer fit of the models. The variation in the 
unadjusted and adjusted data remained unchanged for the 


Richards function. 


At weaning, although the variation in the hean 
prediction error estimates was reduced by adjusting the data 
based on the Richards, Brody and Logistic functions, the 


mean prediction error was only improved in the Richards 


(Table 25). 


At yearling age adjustment of the data reduced the 
variation of the mean prediction errors in ail functions 


significantly improved the fit of the Logistic function 
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Table 26. Standard deviations of mean prediction error by. 
function and age for the unadjusted and adjusted 
data from the Hereford breed group 


SP SP OS OO OW OSS 8 O48 OOS DOB DDS PD PUBOVWlS CHWS TW SBea Qa2 ee qGeeordee ®emwaqoae ooo aoe 


Standard deviations (3%) 


Age Data Richards Brody Bertalanffy Logistic 
ee amaor ee ae ED MY ED GY GQ AD OP aH GP GH) WY ( <P GY GD GP aw OW AP AD GY GD GH a cd GP CY <P aH = <P ap a ah ae a) GE Oe a or Oe we @ @ a OO 
Birth Unadj. 13.184 14.743 17.218 27.2940 
Adj. 13.285 172909 15.346 24.732 
Weaning Unadj. 16.952 © 16.387 14.719 14.662 
Adj. 15.351 15.968 14.529 13.331 
Yearling Unadj. 12.724 12.653 122309 12.265 
Adj. 11.284 11.533 10.922 10.429 
20 Months Unadj. 9.402 9.487 9.807 10.130 
Adj. 10.667 10.505 10. 866 11.337 
32 Months Unadj. 11.888 11.959 12.105 12.252 
Adj. 11.027 11.976 11.037 11.029 
44 Months Unadj. 9.688 9.671 9.501 9.371 
Adj. 9.522 9.595 9.364 9.119 
56 Months Unad4. 8.890 8.801 Be441 8.223 
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(Table 25). 


At 20 months of age adjusting the data increased the 
variation of error in the data while at 32 and 44 months of 
age the variation in the estimate was similar between the 
unadjusted and adjusted data. The variation in the Hereford 
data based on all functions was somewhat reduced by 
adjusting the weights at 56 months but the reduced variation 
was not associated with a reduction in the mean prediction 


errors (Table 25). 


Table 27 shows the mean prediction errors at each age 
based on the Richards, Brody, Bertalanffy and Logistic 


growth functions for the unadjusted and the adjusted data 


from the Hereford population. 


The mean prediction error estimates based on the 
Richards function ranged from -5.9 to 1.7 in the unadjusted 
data, falling into two significance classes while the mean 
prediction errors based on the Brody function showed a 


slightly wider range from -5.9 to 1.9 and fell into three 


Significance classes. The mean prediction error estimates 


based on the Bertalanffy function fell into four 


Significance classes while the Logistic function was least 


consistent and the mean prediction errors classified into 
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four significance classes, with a wider range in the error 


estimates. 


In the adjusted data the Richards means ranged 
from -7.3 to 4.3 falling into four significance classes. 
The mean prediction error estimates based on the Brody 
Showed a still wider range falling into four significance 
classes. Thus, the Richards provided a more consistent fit 
than the Brody across ages. The mean prediction error 
estimates based on the Bertalanffy showed a much wider range 
than one! Richards and Brody falling into five significance 
classes while means based on the Logistic fell into six 


Significance classes. 


qe eee ao aoa 


The Richards was the only model that provided a 
consistently good prediction of birth weight in both sets of 
the Hereford weight data. The Bertalanffy and Logistic 


functions were very unreliable in predicting birth weight 


accurately. 


The Brody function provided a low mean prediction error 
estimate in the unadjusted Hereford data and the Richards 


provided a significantly lower mean prediction error in the 


adjusted data at weaning. 
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Yearling weight was predicted well by all four growth 
functions in the data sets with the Brody function 


estimating it with a consistently low mean prediction error. 


The weight at 20 months of age was estimated accurately 
by the Richards and Brody functions in the unadjusted data 
and in addition by the Bertalanffy function in the adjusted 


data. 


Weight at 32 months was overestimated by all functions 
and none of the growth models could adjust to the temporary 


weight loss. 


At 44 months of age, the Richards, Brody and 
Bertalanffy functions predicted weights with a minimum of 
error and the Bertalanffy provided the lowest mean 


prediction error to the adjusted data. 


At 56 months the Richards and Brody models provided a 
good fit to the unadjusted Hereford data and the mean 
prediction error estimates based on the same functions were 


higher but not significantly different between functions in 


the adjusted data. 


Adjusting the data did not significantly improve the 
fit of the Richards function at birth although the _ nean 
prediction error estimates were reduced. However the fit of 


the Brody function was significantly decreased by adjusting 
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the data and is contrary to what was observed among the 
Hy brids. The fit of the Bertalanffy and Logistic functions 
at birth was improved by adjusting the data. The fit of the 
Richards function was improved by adjusting the data at 
weaning. Adjusting the data had no effect in improving the 
fit of the four growth functions at any of the subsequent 
ages, with the exception of the Bertalanffy and Logistic 


functions at 20 months of age. 


The consistency of fit based on a Duncan's comparison 
of means showed that the Richards function was most 
consistent in estimating mean prediction errors which were 
small and uniform within each set ot Hereford data across 
ages. 


from the _XB_groups 


6.7 Analysis_of data fr ms 


6.7.1 Analysis of the goodness _of fit _of each growth 


function applied to each age 


The mean prediction error estimates, standard 
deviations, standard errors, minima and maxima for the seven 
ages by function for the 100% BR, 75% BR, 50% BR 50% LD and 
75% LB breed groups are shown in Appendix Tables 3 to 6. 
The results show exceptional similarity to the two large 
breed groups with the Richards consistently providing a good 


prediction of birth weight, the Brody providing a less 
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consistent fit at birth and the Bertalanffy and Logistic 
functions providing a poor fit at birth. The weaning 
weights were predicted satisfactorily by the Richards and 
Brody functions while yearling weight was predicted well by 
ail functions in each breed group. Weight at 56 months was 


predicted well by the Richards and Brody functions. 


6.2j/.2 Estimation of the consistency of fit 
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Mean prediction error estimates arranged ina Duncan's 
table derived from the Richards, Brody, Bertalanffy and 
Logistic functions for the 100% BR, 75% BR, 50% BR 50% LD 
and 75% LB breed groups are shown in Appendix Tables 7 to 
10. The Richards function provided a uniform estimation of 
error in the 100% BR, 50% BR 50% LD and 75% LB breed groups 
followed by the Brody which provided an equally consistent 
ET t, In the 75% BR breed group the Brody function followed 
by the Richards provided the most consistent fit across 


ages. 


6.7.3 Conclusions 


The Richards function provided a good estimate of birth 
weight, weaning weight and adult weight while the Brody 
function showed an inconsistent fit at birth and a good fit 


at weaning and 56 months of age. The fit of the four models 
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during the intermediate ages was similar and followed the 
Same patterns as in the larger breed groups. In general, 
the Richards and Brody functions were consistent in 
predicting with small and similar mean prediction error 
estimates across ages which was similar to the observations 


in the two larger breed groups. 


6.8 Correlations between parameter estimates for the 
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Richards_and Brody functions 


The simple corretations between parameters (A), (B), 
(kK), (4), (AGR6), (AGR12), (AGR18) and (¥*) from the 
Richards model fitted to the Hybrid and Hereford lines are 
Shown in Table 28 and for the Brody model in Table 29. 


Correlations were computed using adjusted data. 


Correlations between (A) and (k) based on the Richards 
and Brody growth functions were similar in both Hybrid and 
Hereford breed groups ranging from -0.53 to -0.66 indicating 
early maturing animals grew to smaller adult weights. 
Similar relationships were reported by Brown et al. (1976) 
in Jersey and Hereford females, and by Fitzhugh and Taylor 
(1971) with Herefords. Thus, selection for early maturing 
types should lead to correlated responses in reducing adult 
body weight. In the present analysis the maximum difference 
in correlations between parameters (A) and (k) in the two 


breed groups based on the Richards and Brody functions was 


Pe vide trawo les bas ak 
vs "cheno be ae 
‘me, cf wad 


e | aa its dw | 


w 
if 


= 


er way ate itn ey envi sein A 

si oh ee Bas orb } 

wie band nets | “ie Me. ‘of bared 
ies: ered ae tahoa hod ee oy | bas oo: i 


- 
i 


on ‘at: ec, ua ata 


seein, Dae pede 
veedgnen. 2 inbe. aot haie S 
(avery ste £5. hain hae 
WiC oR “pga! ato |e 
pmukatiew a) Ed was i ai 
ales eiidintin snd ic 

~ epndaeta.tp ere | 

7 1 anid ay brie ‘nilak | 


[= 


r 


LO°O>d ex 
S0°0>d* 
= Cl o= *xES°O xen °O *k08°O **O0L °O eciine Lz*0- Ged) 
#Bh °O = Peryar wet °O xB E°O *x OE °O x66 °O= LA) (gLuoV) 
*nl °O men B ° 0 = xxhB°O #x9L°O **08°O **9L °O~ Ec 0= (ZLuoW) 
*€9°O xx 19 °O x76 °O = *xlG°O *%*x9L °O #x 9S °O= LL°0= (9u9¥) 
x16 °O aLE°O | **BL°O *%08°O - eeiB°O  #**66°0= eS 0- (Ww) 
*Z7S°O *9L°O **99°O *xx08°O xe b6°O - #x BB °O= #*#x99°0= (x) 
*& 16 °O= *¥LE°O= **LL°O= *xEL Om *xB 6° O= *%98 °O= - xa Lt ?0 (a) 
9€°0 walt °O f0°O c1.70= xaGn°’ Qe #*xG9° De] *xbLp°o * (vy) 


SOS SS OS OSSSSSSSSlFQeLeS SRBSSBZ®L®®2®® 2S eheowre> Se BSSe @F SSOP SSS SSS OF SSSSSSeS@MOOeest © VPeE@eeeBeeena aoe ee 


(xX) igLugv) (ziyow) (ou9y) (nw) (4) (a) (v) 


SSS SEG E688 SSO SOSSE SSHeSOO2Seovt eV SSPVST SST SSBO® See te SETS eSe Goes Ve Ser Ceo SSO SSP eee ee © V2SGeese’ 2e22Teee 


(yeuobetp aoteq) ejyep peysnfpe prozyezey pue (Tewobetp aaoqe) eyep peysnfpe 
praqky ay2 0F pez3TF Tepow spreyotyY e4y worZ Slazemered yyAoIHh uwssazeq SuOTJeTeTIOD °szZ ETqeL 


- : 
vy : 
AP 
iz 
mh ms 
Py, 4] 4) 7 
‘ 
ri ’ 
ims 
ail 
rs 
i 
De 
A : 
0 
’ 
t 
A 
t 
3 
} 
pl, 
7 
7 
i 
Md 
t 
ve 
ie eS 
= y 
e 7 1 
* 
i 


118 


LO°O>d xx 

G0°O0>dx 

OS OoT Se See Qe GS SSG ee SE OOo se eer eo See ese eee SF eoeoeeeewde @®& @& SFeeese @ee eo @ooe SS Oe e2ee oe eee @ eee ce & @& 
= ae EB °O 6L°O #xEL POS 70°0 *x96°O (gLuowv) 

xe16°O ~ ##O2°O * EZ PO ZL°O “xx 8B°O (ZLuoV) 
eLicf°O wxxZ i °O = ee ES °O 9L°O x hE °O {ou9V) 

xeEL °O= *67°O= *xGh°O = 40°O xx BGP O= (») 
60°O *xBZ°O *07°O 6L°O - EL°Oe (a) 

ex S6°O ex EO 6 °O eB oO ex EG Oe 9£°0 ~ (yw) 


R2PSSOS SSS ® SO&SS®eSGoe a hae on @OeeoS SS OSS S28 S288 SSS S228 S65 8 SSS SSES2SSSSSE EUR OOUEee 


(gLyoyv) (ZL aowv) (ox9¥) (x) (a) (y) 


SBSSSD 2 SSeS VS SSWSBWS*e SF SSO SEOoo SOG 6 SBSSES SSS STS BSSSSOSSOSSSSVt>OQ eSesecoo eon @oe 


(Teuobetp aojtTeq) eirep proyezey pe ysn€pe pue 
(Teuobetp eaoqe) erep ptaqdky peasn€pe ayy 0% paz3TJ Taepou 
Aporlg 343 worzy sTeZouerted yIAoIh usaanzyaq SuOTZeTAeIIODN °6Z2 STaeL 


<i 


¥ 


v r | 
A 
eas 


119 


0.05. In a preliminary analysis, using the unadjusted 
weight data, the difference in the correlations between the 
parameters (A) and (k) based on the Richards function 
between lines was 0.14 (see Appendix Table th 

Furthermore, the % difference between the maturing rates (k) 
of the Hybrid and Hereford lines based on the Richards 
function was 17.1% for the adjusted weights (Table 18) while 
there was an 8.82% difference in their asymptotes. Brown et 
al (1976) observed a) wider differences between correlations 
of (A) and (k) between breeds and b) a greater percentage 
difference in the maturing rates compared to the asymptotes 
between two breeds and they were of the opinion that the 
parameters (A) and (k) are relatively independent. However 
from the results of the present analysis there is little 


evidence to suggest that the two parameters are independent. 


In the Hybrid breed group correlations between (A) and 
{(AGR) based on the Richards function changed from negative 
to positive (r=0.47) as the age increased from 6 to 18 
months. A similar trend was observed in the Hereford breed 
group although, the correlation of (A) and (AGR18) was 
Slightly less positive compared to the Hybrids. Thus, 
animals with lower than average growth rates during the post 
weaning to yearling ages (6 to 12 months) grew faster 
relative to its mature weight, and animals with higher 


growth rates at 18 months were again heavier at maturity. 
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Correlations between (A) and (AGR) based on the Brody 
function fitted to the two breeds were Significant and 
positive at 6, 12 and 18 months of age, with the values 
increasing with age. Animals with higher growth rates at 
later ages were therefore more likely to achieve higher 
mature weights than animals with higher growth rates at 


early ages. 


Correlations between (A) and (AGR6) based on the 
Richards function were negative (r = -0.12 and -0.17) for 
the adjusted Hybrid and Hereford data, respectively, 
although the values were not significant. However, in both 
breed groups (AGR6) was positively correlated with the 
maturing rate parameter (k) (r=0.80 and r=0.76) suggesting 
that early maturing animals with higher gains at 6 months of 
age grow to smaller mature weights. The same effect was 
less pronounced in the Hybrid line when gain at 12 months of 
age was considered, as the correlations between (A) and 
(AGR12) were less negative and the correlations between (k) 
and {AGR12) less positive. The relationships were reversed 
at 18 months of age as the correlations between (A) and 
(AGR18) were positive (r=0.47 and r=Q.17) in “Hystidy and 
Hereford lines, while the correlations between (AGR18) and 
(k) were less positive (r=0.16 and r=0.30) in the same two 
lines. The correlations between (A), (AGR18) and = (k) 


suggested that late maturing types with higher growth rates 
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at 18 months of age had higher mature or adult weights. 
Brown et al (1976) found early maturing types with higher 
growth rates at the point of inflection grew to smaller 


Mature weights. 


Based on the Brody function, animals with a slower 
general maturing rate and rapid growth rates after 1 year of 


age were heavier at maturity. 


According to the formulae of the models when (M) 
increases (Y¥*) should increase and this is shown by the high 
and positive correlations between the two parameters. Pe a 
would therefore suggest that populations which are 
Characterised by higher absolute values for (M) are heavier 


at inflection. 


The correlations between (A) and (M) were negative with 
values of “0.45 in both breed groups. Thus, the 
correlations between (Y*) and (A) was variabie between 
breeds, ranging from -0.21 in the Hereford breed to 0.36 in 
the Hybrid breed. A positive correlation between (Y¥*) and 
(A) in the Hybrid line suggested that heavier animals at 
inflection were heavier at maturity. However, the 
correlation between (Y*) and (A) was not significantly 
different from zero. The correlations between (Y*) and (A) 
May again be influenced by the absolute values of (A), as 


(Y*) is a function of (A) and (M). Fitzhugh (1976) stated 
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that it is often difficult to interpret relationships based 
on algebraic models as characters are not completely 
independent, and often relationships contrary to a basic 


biological trend are obtained. 


The correlations between (k) and (AGR6) were high, 
ranging from 0.76-0.89 in both breed groups based on the 
Richards and Brody functions, suggesting that early maturing 
types grow rapidly to 6 months of age. The correlations 
between (k) and (AGR12) based on the Richards function were 
less positive in both breed groups. The correlations 
between (k) and (AGR12) based on the Brody function and 
applied to the two breed groups was negative. Based on the 
Richards function, in the Hybrid and Hereford lines the 
correlations between (k) and (AGR18) were slightly positive 
with values of 0.17 and 0.30 respectively. At 18 months the 
correlations between (k) and (AGR) based on the Brody 
function applied to the breeds were negative indicating that 
at 18 months early maturing types had lower absolute daily 
gains. In general, early maturing animals were 
characterised by faster growth rates during early ages and 


slower growth rates at later ages. 


Weight at the point of inflection (Y*) was positively 
correlated with (AGR6) in all of the data and negatively 
correlated with (AGR18) in the Hereford line suggesting that 


Herefords heavier at inflection had higher growth rates at 6 
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months of age but lower growth rates at 18 months of age. 


The same result was not observed among the Hybrids. 


The correlations between (M) and (k) were high and 
positive along with the correlations between (M) and (Y*). 
This would suggest that early maturity was thus associated 


with higher values for (M) and higher inflection weights. 


The correlations between (M) and (AGR) at 6, 12 and 18 
months showed the same trend as for the correlations between 
(kK) and (AGR) within each breed group. Thus, heavier 
animals at inflection showed higher growth rates at 6 months 
of age but relatively lower growth rates at 18 months of 


age. 


Interrelationships between growth rates at different 
ages have often been investigated on actual data and 
Fitzhugh and Taylor (1971) observed that a higher growth 
rate during one period was accompanied by a lower growth 
rate in another. Based on the Richards model, in both the 
Hybrid and Hereford data the correlations between (AGR6) and 
(AGR12) were high, ranging from 0.94 to 0.79. The 
correlations between (AGR6) and {AGR12) based on the Brody 
model were lower ranging from 0.72 to 0.39. Similarly the 
correlations between (AGR12) and (AGR18) were always 
positive; the correlations based on the Brody model being 


somewhat higher. Thus, selection for growth rate at 6 
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months should increase growth rate at 12 months, and 
selection for an increased growth rate at 12 months should 
increase growth rate at 18 months of age. However, as the 
age interval increases the degree of association decreased. 

Thus, when ages are widely separated, selection for a high 
growth rate at the lower age would not be expected to 


rapidly improve growth rate at the higher age. 


6.8.21 Conclusion 


In all derived biological models the relationship 
between the maturing rate and mature weight was negative 
suggesting that lighter animals at maturity are early 
maturing. Barly Maturing types had absolute growth rates 
that were above average at early ages (6 months) and below 
average at later ages (18 months), suggesting that selection 


for growth rate at early ages would increase maturing rate. 


Based on correlations of the Richards model animals 
with higher adult body weights grew more slowly at 6 and 12 
months of age and grew faster than average at 18 months 
although growth rates were relative to adult body weight. 
Based on the Brody model, animals with higher growth rates 
after 12 months were expected to be heavier at maturity. 
Furthermore, early maturing animals with higher gains at 6 


months grew to smaller mature weights and late maturing 
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animals with higher gains at 18 months grew to larger mature 


weights. 


Absolute growth rates were positively correlated when 
the two ages were closely adjacent but less positively 
correlated when widely separated suggesting that selection 
for increased growth rates would not necessarily improve 
growth rate to the same extent at subsequent ages especially 


if the two ages are widely separated, 


6.9 Analysis of absolute growth rates 


Table 30 contains the observed absolute growth rates 
(AGR), the predicted absolute growth rates based on each 
growth function, and the differences between the observed 
and predicted growth rates by growth period and breed group 


for the adjusted data from the Hybrid and Hereford lines. 


Figures 9 and 10 show the observed and predicted 
absolute growth rates by function and breed group for the 


adjusted data. 


The highest observed absolute growth rates (average 
daily gain) were recorded during the pre-weaning stage in 
both breed groups; post-weaning to yearling gains were 
somewhat lower than the pre-weaning gains and the drop was 
greater among the Hybrid heifers. The growth rates 


increased again during the yearling to 18 month period in 
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both groups and dropped to a minimun during the stage from 
18 months to adult. The cascading effect seen to occur in 
growth rate between adjacent periods was also observed by 
Fitzhugh and Taylor (1971). The absolute growth rates 
predicted by the Richards, Brody, Bertalanffy and Logistic 
functions were significantly different from the observed 
(P<0.01). In both breed groups the variations in season and 
Management accounting for weight changes limited the use of 


growth functions in predicting absolute growth rates. 


During the period from birth to weaning the Richards 
function predicted mean absolute growth rates with a 4.7% 
and 2.2% error in relation to the observed growth rates in 
the Hybrid and Hereford lines respectively while the Brody 
function predicted absolute growth rates with a 4.1% and 
6.6% error in the same lines. The Bertalanffy and Logistic 
functions were poor in predicting gains during the pre- 
weaning stage due to an overestimation of birth weight and 
an underestimation of weaning weight. All growth functions, 
with the exception of the Brody appiied to the Hereford 
line, underestimated pre-weaning absolute gain. The 
overestimation based on the Brody function of gain in the 
Hereford line was due to an underestimation of birth weight 
and an overestimation of weaning weight which increased the 
difference between the observed and predicted weights. In 


the Hereford line the closest approximation of the mean 
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observed pre-weaning gain was seen with the Richards 
function when birth weight was slightly underestimated and 


weaning weight was slightly overestimated. 


During the post-weaning to yearling period all 
functions applied to the Hybrid data overestimated absolute 
growth rates; the Richards by. (20.6%, Brody) 202%, 
Bertalanffy 37.2% and the Logistic by 36.2%. The growth 
rate predicted by the Brody function for the Hereford data 
fitted very well. Furthermore, as there appeared a general 
overestimation of daily gains during this period, the growth 
rates computed by each of the functions were closer to the 
observed growth rates among the Hereford heifers as they 
showed higher growth rates during the feed test period which 


followed weaning. 


The deviations of the predicted gains from the observed 
were high during the yearling to 18 month period in both 
lines; observed daily gains were also higher compared to the 
preceeding stage in both lines as the period coincided with 
summer when the heifers generally showed greater increases 
in body weight. Due to this temporary environmental effect 
that increased weight gain over the period, greater 
differences were observed between the predicted and observed 


gains, all functions underestimating growth rates in both 


breeds. 
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The predicted daily gains during the last stage were in 
closer agreement with the observed daily gains as growth 


rate was more uniform and the time interval large. 


A similar analysis was done for mean gain over period 
comparing the observed weight gain and the expected weight 
gain using each growth function. Appendix Table 12 shows 
the mean weight gain over the period and the differences 
between observed and predicted weight gains by breed group 
and function. The results were similar to the analysis of 


absolute growth rates. 


As growth functions regress size on age and are of an 
exponential decay type, they are poor in adjusting for 
fluctuations conditioned by Management practices and 
environment, often seen in biological data. The gains 
predicted on these models assume a progressive decline in 
growth rate at each period and thus the two functions that 
provided a good fit to the data, i.e. Richards and Brody 
almost follow a straight line across all periods (Figures 9 
and 10). A display of the observed gains by a linear 
regression of growth rate on age with the seasonal effects 
removed would have afforded a worthwhile but more 
theoretical appraisal reflecting the genetic trends in 
absolute growth rates of the population. Furthermore, by 


this technique the observed gains at each period could have 
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been more appropriately compared with the predicted growth 


rates based on the 4 growth models. 


The Richards and Brody functions provided estimates of 
predicted preweaning absolute growth rates to the adjusted 
data from the Hybrid and Hereford lines which limited the 
mean error to <7%. The Brody function was particularly 
efficient in predicting mean absolute gains similar to the 
gains observed in the Hereford line during the weaning to 
yearling period. All growth models were poor in predicting 
absolute gains during the yearling to 18 month period. The 
absolute gains predicted by the four growth functions were 
quite similar to the observed gains during the last growth 
period even though, only the predicted absolute gains based 
on the Bertalanffy function were not significantly different 


from the observed absolute gains in both breed groups. 
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7e OVERVIEW 


The four parameter Richards model as well as the 
Simpler 3 parameter Brody model fitted all of the data 
better than the Bertalanffy and Logistic models. The 
Bertalanffy and Logistic models provided a poor fit as birth 
weight was always overestimated and adult weight 
underestimated. However, none of the functions could 


account for seasonal fluctuations that reduced body weight. 


A growth curve that provides a good fit toa set of 
longitudinal observations of a single animal or population 
reflects the genotype of the animal or population and 
selection for growth traits on a predicted growth curve may 
be carried out if the curve is adjusted for temporary 
environmental conditions. In any animal improvement program 
which is founded on the principles of breeding and selection 
as a basis of improvement, selection for growth traits 
reflect changes in the shape of the basic growth curve. 
Selection for pre and post weaning growth rate, which is a 
common practice in essence selects positively for an 
increased maturing rate. Furthermore, growth rate and body 
weight within a stage of growth is positively correlated 
especially at early ages (Fitzhugh and -Taylor, 1971). By 
selecting for increased growth rates prior to 1 year of age, 
weight at that age and maturing rate are increased while 


mature weight may be decreased. However the increased 
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growth rate is only relative to mature size. If the effects 
of selecting for increased growth rates were to be shown on 
a cumulative growth curve, the slope of the curve prior to 
inflection would be increased and the Slope post inflection 
decreased. From this discussion it is very evident that the 
time at which beef animals are marketed becomes important, 
in order to take advantage of these growth characteristics, 
as it becomes less economical to keep animals well past 


their peak growth periods. 


In the recent past the effect of cow size and 
production efficiency has been emphasized (Long et al., 
1975a; Morris and Wilton, 1976). Selection for an increased 
Maturing rate may be one method of reducing size as it is 
negatively correlated with mature weight. However, as_ the 
Mature body size or adult weight is highly heritable unless 
selection is intense the rate of progress would be slow. 
Theoretically, selection for increased adult body size on 
the other hand will have opposite effects and would select 
late maturing types in which the change of the slope of the 
cumulative growth curve would be more gradual. It must be, 
however, born in mind that the preceeding discussion is 
limited to selection within a breed and should not be 


confused with slow or fast maturing breeds of cattle. 


Indirect selection for increased maturing rates can 


also have an effect on carcass quality. Ideally, a _ beef 
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animal at the time of marketing should contain maxinun 
muscle, and optimum fat and minimum bone (Berg and 
Butterfield, 1976). Animals at a higher degree of maturity 
would be more likely to put on fat compared to muscle, as 
muscle growth reaches a plateau at a smaller weight and age 
while fat growth can increase until the animal approaches 
its adult body size at which time it will plateau. Once the 
animal reaches adult body size which is under genetic 


control, the relative gain will be small. 


A series of weight-age points of a longitudinal type is 
often difficult to interpret; one method of condensing all 
of the information into biologically interpretable 
parameters is by the use of growth curves. When the values 
of the fitted parameters are known, many of the commonly 
used growth traits such as average daily gains, 
instantaneous growth rates, lifetime growth rates, relative 
growth rates and maturing rates could be easily obtained. 
However, the accuracy of the defives estimates such as 
growth rates using the growth ‘constants (Brody, 1945) or 
fitted parameters depends on how well the function of choice 
fits the data at each of the selected data points. Once the 
function of choice has been identified the data generated by 
the function can be used as an alternative to the true 


biological model as if it was known (Brown et alse, (12:26) 6 


An appropriate method of estimating the fit of a growth 
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model is still under discussion. Fitzhugh... (1976) states 
that the methods such as the use of residual variances may 
be subject to error due to correlations between longitudinal 
data points. Furthermore, the correlated effects may be 
stronger if seasonal effects are significant in determining 
weight at certain times. Allen (1967) as reported by 
Fitzhugh (1976) described a nuiltivariate procedure 
appropriate for non-linear models assuming that observations 
on different animals at the same age were correlated. 

However, in data of the type analysed in this study such 
assumptions may not hold and in the analysis of biological 
data removing the correlated errors between longitudinal 


data points is difficult. 


In the present analysis a pooled residual variance and 
mean prediction error estimate were used to calculate the 
fit of the functions assuming that the correlated errors 
between ages were removed in the weights predicted by the 
growth functions. The usual test of fit such as the sum of 
Squares of the deviations (y=) 2 or a mean square deviation 
was not used due to the inherent bias existing in this 
method of analysis which attributed an equal weighting when 
the observed weight was 50 kg and predicted weight 100 kg or 
when the observed weight was 500 kg and predicted weight 550 
kg. Thus, a procedure using a relative fit or mean 


prediction error estimate was considered appropriate when 


mene ahha, Eas: tite 


oth ven asdialxie | os isin Le | | | ae AF 
adh bay Says ded) ak vawonanbes enaiesonie. | ine 


ee Peete agers te RE “tht 2oubry ‘a ; 7 ion ” 
zy bi Een | eons i” \begraoe elas serene i Anas 
ooh tes nee ae tga eh Wawen | ell | wen tAonae ponies ad P 

Bay 5 26 1 tere +), “ age eee ae | “ig ehentng 1h ' ” 
dune Outs Fa aes ih Septet Git ar ‘to ey Go 


Caimitye Toke sor iets withs ‘S48 ae tds S104 ohm 


; 


De. al 


‘Ce ef ny 
f wi stipe. nkaerod, ans oats Le inapea ear ie : 
; ek ' if a ¥ ae odthe 
7 Ane. ha he ie ree : 


adit shaun sion tony wae aoe ys 


naoute pavsiteas OD oehe aig lnendws sania tt 
aid yet besaiiaag pas hiy ava a fy apna Bi! meee. 
to’, wpe, esis as 426R nit Oe va Kee aa 


Pibtgh Be Hh Pe or 


Ho Ltk owed lala Se fie #: ‘eee * am 


f Kd. +e AW, 


a a a Neola 
mrad ag nse) ‘itl Soe at | ‘wae: cl bass: ne 
. : ih as : ha mei, ae geet 
A serine bande | body 
tio ge ony, feubow sassesonl fie’ | 
ew ‘iy i 


jee paiphay ‘pale Beg 3 
bigs: +4, atee Beis # 


CcPashi paths ages sit pnp 


137 


the fit of the four functions was to be evaluated. Brown et 
al. (1976) and Fitzhugh (1976) in addition to using 
residual variances, relied on visual comparisons of tae 
fitted curves with the actual data at selected ages together 
with the biological interpretability of the fitted 


parameters aS a measure of the accuracy of fit. 


In studies comparing fitted parameters derived fron 
growth functions applied to two breeds or sets of data the 
importance of the ccnsistency in the error estimate has not 
been stressed in the literature. Thus, in an ideal 
Situation not only must the function fit the two sets of 
data well but errors calculated at each age or overall ages 
must be consistent in order to make a reliable comparison of 
the fitted parameters such as (A) and (k) especially as they 
are negatively correlated. The present analysis showed the 
importance of the consistency of error when estimates of the 
same parameter based on a biological model is compared 


between two breeds or sets of data. 


A knowledge of growth curves is important to scientists 
interested in lifetime trends. The application of growth 
curves to a part of the lifetime data does not provide 
accurate information about the growth patterns of an animal 
or population. Weight age data on each animal from birth to 


an appropriate adult weight should be available for a growth 


curve analysis. 
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In conclusion growth models offer a Mathematical 
description of growth. The generalized Richards or the 
Simpler Brody modei can permit a detailed study of growth 


and its properties in beef cattle. 
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Appendix Table 11. Correlations between growth 
parameters from the Richards 
model fitted to the unadjusted 
Hybrid data (above diagonal) 
and Hereford data (below 


diagonal) 
(A) (B) (k) (M1) (AGR12) 
(A) — 0.51 -0./70 -0.51 “0.15 
(B) 0.56 faa -0.85 -0.97 -0.56 
(kK) -0.86 -0./78 ru 0.89 0.66 
(M) -0.62 -0.97 0.84 ax 0.68 
(AGR12) -Q.44 -0.56 0.7 0.67 eed 


All correlations were significantly different 
P<0.01 
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